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ABSTRACT 

The  experimental  and  theoretical  results  of  an  optically  pumped  vibrational  CO 
laser  are  explored.  A  tunable  1550-1580nm  optical  parametric  oscillator  laser  was  used 
as  a  pulsed  pump  source.  The  OPO  pumped  the  R(6)  through  R(9)  (3,0)  overtone  band 
of  the  CO,  which  induced  lasing  on  the  (3,2)  and  (2,1)  bands  around  4.7um.  The  laser 
output  was  spectrally  separated  to  determine  the  spectral  and  temporal  evolution  of  the 
CO  lasing  pulse.  The  influence  of  CO  pressure,  cell  length,  pump  pulse  intensity,  and  the 
resonator  Q-Factor  on  the  CO  laser  pulse  temporal  profile  and  spectral  distribution  were 
studied.  A  numerical  model  was  developed,  and  the  experimental  results  are  compared 
to  the  model  predictions. 
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CHAPTER  ONE:  INTRODUCTION 


1.1  The  Optically  Pumped  Carbon  Monoxide  Laser 

This  dissertation  describes  the  most  extensive  investigation  into  a  2nd  overtone 
optically  pumped  Carbon  Monoxide,  CO,  laser  to  date.  The  2nd  overtone  CO  laser 
consists  of  a  pulsed  pump  source  operating  at  approximately  1.5  fan  used  to  exicte 
various  rotational  transitions  on  the  2nd  overtone  band  of  CO.  The  CO  lases  twice  on  two 
successive  inversions,  first  from  the  (3,2)  band  and  then  from  the  (2,1)  band  both  around 
4.7  fan .  This  thesis  is  intended  to  adequately  describe  the  radiative  and  kinetic  processes 
that  affect  the  optically  pumped  CO  laser  and  to  shed  light  on  future  improvements  to  the 
4  micron  CO  cascade  laser. 

Carbon  monoxide  has  been  developed  as  a  high  power  electric  discharge  laser;  it 
was  first  developed  in  1964  at  Bell  Laboratories  [1].  The  electric-discharge  CO  lasers 
operate  on  a  many  high  lying  fundamental  bands,  and  exhibit  broad  band  output  around 
5-6pm  [2]-[3].  However,  the  atmospheric  absorption  window  is  closed  in  the  output 
region  of  the  electric  discharge  CO  laser,  see  Figure  1.1.  Optically  pumping  CO  provides 
a  mean  to  excite  a  specific  rovibrational  level,  and  lase  on  a  single  fundamental  band  at  a 
wavelength  around  4.7pm,  which  falls  within  the  atmospheric  window.  Another  useful 
application  of  optically  pumping  the  2nd  overtone  of  CO  is  the  convenient  wavelength 
required,  1560- 1576pm.  This  wavelength  happens  to  coincide  with  Erbium  fiber  laser 
wavelengths,  which  have  been  developed  for  the  telecommunication  industry.  Optically 
pumping  CO  shows  potential  as  an  efficient  means  to  beam  combine  and  frequency 
convert  the  output  of  a  high  power  fiber  laser  bundle  to  a  single  4.76pm  beam  with  good 
beam  quality. 
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Various  aspects  of  the  CO  lasing  process  were  studied  experimentally  and 
theoretically  and  will  be  presented  in  this  thesis.  The  influence  of  CO  pressure,  cell 
length,  pump  pulse  intensity,  and  the  resonator  Q-Factor  on  the  CO  laser  pulse  temporal 
profile  and  spectral  distribution  were  studied.  A  numerical  model  was  developed,  and 
the  experimental  results  are  compared  to  the  model  predictions.  At  longer  gain  and 
resonator  lengths  (2m)  the  CO  laser  ran  in  mode  locked  operation. 


1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0  5.5  6.0 

Wavelength  (|jm) 


Figure  1.1.  Atmospheric  transmission  spectrum.  Notice  the  high  transmission  around 
4.6-4. 8um. 

Optically  pumping  a  molecular  gas  shows  potential  in  becoming  a  viable  means  to 
produce  a  high  power  4.6  /Jm  laser  system.  The  wavelength  region  around  4.6-4. 8  /um  is 
of  great  interest  due  to  the  atmospheric  window  in  that  region,  Figure  1.1. 

1.2  Previous  Studies 

Carbon  monoxide  has  been  used  as  an  optically  pumped  lasing  medium  in  the 
past.  The  experiments  conducted  at  AFRL  demonstrated  CO’s  ability  to  be  optically 
pumped  at  2.36pm  and  lase  from  the  1st  overtone  band  [4]-[5].  McCord  demonstrated 
multi-line  lasing  in  (2,1)  band  of  CO  with  a  reported  efficiency  of  25%  in  1999  [4], 
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Further  work  by  McCord  showed  single  line  lasing,  with  a  reported  efficiency  of  13% 
from  the  (2,1)  band  by  introducing  a  grating  into  the  resonator  [5]. 

Various  other  molecular  gasses  such  as  HBr  and  HC1,  were  optically  pumped  to 
higher  vibrational  overtones  [6]-[8].  Both  gasses  successfully  produced  cascade  laser 
outputs  around  4  pm.  The  above  mentioned  research  was  referenced  heavily  for 
experimental  approaches  and  considerations. 

1.3  Thesis  Layout 

Chapter  2  reviews  the  various  theoretical  concepts  required  to  adequately 
understand  and  model  this  CO  laser.  A  brief  discussion  of  the  following  topics  are 
included  to  refresh  the  reader;  diatomic  vibration  rotational  energy,  kinetic  processes, 
radiative  processes,  transverse  and  axial  resonator  modes,  and  mode  locking.  The 
discussions  are  not  intended  to  be  complete  in  their  theoretical  derivation,  but  rather  act 
as  a  basic  guide  to  the  physics  used  in  the  data  analysis. 

Chapter  3  outlines  the  experimental  setup  and  procedures  used  in  the  laser 
experiments  in  great  detail.  All  the  instruments  used  in  the  data  collection  will  be 
adequately  documented.  The  operating  and  aligning  procedures  for  the  optical 
parametric  oscillator  are  included  as  a  reference  for  future  research. 

The  results  and  conclusions  of  both  the  theoretical  and  experimental  data  will  be 
presented  in  Chapter  4.  Chapter  5  contains  a  brief  discussion  of  a  future  research  effort 
in  continuous  wave  pumping  of  CO. 
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CHAPTER  TWO:  THEORY  OF  THE  CO  CASCADE  LASER 


2.1  Rotational  /  Vibrational  Energy  Levels 

The  energy  levels  of  a  quantum  mechanical  system,  are  defined  by  the  solution  to 
the  Schrodinger  wave  equation. 


=  emvm 


(2.1) 


Where  H  is  the  Hamiltonian  operator,  ^represents  a  molecule’s  wave  function  and 
Em  is  the  total  energy  of  the  system  including  electronic,  vibrational,  rotational,  and 
translational  energy  [9],  The  Hamiltonian  of  a  hydrogen  molecule  is: 


"=^77(V>V»)-^1(V?+V5)- 


h 2 


2  M 


+  - 


2  m„ 


e2 

'  1 

4  ne0 

_rM 

•  +  —  +  —  +  — 
r,B  r,,  r. 


4  7TSn 


J_  1 

r\2  E 


(2.2) 


where  M  is  the  mass  of  the  hydrogen  nucleus,  me  is  the  mass  of  an  electron,  subscripts  1 
and  2  refer  to  the  two  electrons,  subscripts  A  and  B  refer  to  the  two  hydrogen  nuclei,  V 
si  the  Laplacian  operator,  and  the  distances  for  rtJ  are  found  in  Figure  2.1. 


Figure  2.1.  Distances  used  in  the  hydrogen  molecule  Hamiltonian. 
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The  first  four  terms  represent  the  kinetic  energy  of  both  nuclei  and  electrons  respectively. 
The  next  four  terms  represent  the  potential  energy  from  electron-nuclei  interactions.  The 
second  to  last  term  is  the  repulsion  potential  energy  between  the  two  electrons,  and  the 
last  term  is  the  repulsion  potential  energy  between  the  nuclei. 

The  hydrogen  molecule  is  used  as  an  example  for  its  simplicity,  and  the  more 
general  Hamiltonian  for  any  diatomic  molecule  can  be  easily  developed  using  the 
hydrogen  molecule  as  an  example. 


H= 


■ h 2 


-r 


v2.  -- ^-V2  -  — -Z— +  (2-3) 


2M  j  A  2  M/b  2me*- 


Z  nC  £ 

•  + 


2  zAzB 


i  rAi  i  rBi  U  riJ 


R 


The  summations  in  equation  2.3  allow  for  an  indefinite  amount  of  electrons  in  the 
diatomic  molecule,  and  Z  is  the  nucleus  charge. 

In  order  to  solve  the  differential  equation,  a  few  substitutions  and  assumptions 
must  be  made.  First,  the  wave  equation  is  separated  into  two  parts,  nuclear  and 
electronic: 

M)  (2-4> 


where  'Fw(/?)  is  the  nuclear  portion,  and  {r,R)  is  the  electronic  portion.  Substituting 
equation  2.4  into  2.1,  applying  the  common  vector  calculus  identity, 

V2  [fg]  =  fV2g  +  2V/Vg  +  gV2f  (2.5) 

and  using  the  Bom  Oppenheimer  approximation  drastically  simplifies  the  Schrodinger 
equation.  The  Bom  Oppenheimer  approximation  simply  states  that  electrons  respond 
instantaneously  with  respect  to  nuclear  motion  since  the  mass  of  an  electron  is  much  less 
than  that  of  a  nucleus  [10].  Also,  to  first  order  the  electron  forces  have  no  effect  on  the 
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nuclear  energy.  This  approximation  forces  all  differential  terms  involving  both  the 
electron  and  nuclear  energies  to  have  a  value  of  zero,  and  can  be  written  as: 


V  -  0  VvF=0 

V  A,B T  el  U  V  i  *N  U 


v2  'F  =  0 

v  A,B  *  el  u 


vX=o 


(2.6) 


The  Schrodinger  equation  can  now  be  written  as: 


+[»'„„+r„+F,Jv K 


The  equation  is  now  separable  into  nuclear  and  electronic  parts. 


(2.7) 


2we  ■ 


IVfT„+(F„+F„+F„,)  =  CT,; 


(2.8) 


2 


1  ,  1  , 

— v2vf„+— v2^ 
ma  a  n  mb  b  n 


=  -CVn  +  Em'¥n 


The  separation  constant,  C,  represents  the  total  electronic  energy  in  a  diatomic  molecule, 
and  is  a  function  of  the  intemuclear  separation,  U(r).  The  total  nuclear  energy  of  the 
molecule  can  be  solved  by  further  separating  the  wave  function  into  angular  and  radial 
parts:  xi’{R,6,(f)=  R(r)Y(d,</)).  The  solutions  to  the  angular  part  of  the  nuclear 
Schrodinger  equation  are  the  rigid  rotor  wave  functions,  which  are  called  spherical 
harmonics  [11].  The  rotational  energy  levels  of  a  diatomic  molecule  using  the  rigid  rotor 
model  are  given  by: 


Ej  =~2jJ(J  +  1)  ^  -0,1,  2,  ... 


(2.9) 
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where  I  is  the  molecule’s  moment  of  inertia  using  the  center  of  mass  approximation. 
Each  energy  level  has  a  degeneracy  of  gj  =2.7  +  1 ,  and  only  AJ  =  +1  transitions  are 
allowed.  The  selection  rules  for  rotational  transitions  will  be  explored  at  a  later  time. 

The  radial  portion,  R(r),  of  the  nuclear  Schrodinger  equation  results  in  vibrational 
energy  levels,  and  the  form  of  U(r)  must  be  known.  Typically,  U(r)  can  be  represented 
by  a  parabolic  oscillator  for  a  first  order  approximation.  The  Taylor  series  expansion 
around  r<j,  where  ro  is  the  equilibrium  bond  length,  is: 

U^)  =  U(h)+U\rJr-r,)+^U\rJr-rJ+^-(r,\r-rJ  +...  (2.10) 

The  potential  energy  is  zero,t/(r0)  =  0  at  the  equilibrium  bond  length.  By  definition  the 
first  derivative  at  r0  is  also  equal  to  zero,  since  the  potential  energy  curve  is  at  a  minimum 
at  ro.  Using  the  first  order  approximation  for  the  potential  energy  function,  one  can  easily 
see  it  is  a  parabolic  harmonic  oscillator  function,  see  Figure  2.2. 

CO  Potential  Energy  Harmonic  8.  Anharmonic 


Figure  2.2.  Potential  well,  harmonic  &  anharmonic. 
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The  solution  to  the  radial  portion  of  the  nuclear  Schrodinger  equation  results  in  Hermite 
polynomials  and  discrete  vibrational  energy  levels. 

Ev=ae{v  +  yi)  v  =  0, 1,  2, ...  (2.11) 

coe  is  the  zero  point  energy,  and  vibrational  transitions  are  limited  to  Av  =  ±l  by 
selections  rules. 

Spectroscopy  Notation 

A  brief  discussion  of  the  spectroscopy  notation  used  through  out  this  thesis  is  in 
order  here.  A  list  of  CO  spectroscopic  constants  can  be  found  in  Table  2.1  published  by 
Dunham  [12],  The  importance  of  each  term  will  be  developed  in  the  following  sections. 
A  P-branch  transition  is  defined  as  A J  =  JUpper  -JLower  =  -1 ,  and  an  R-branch  transition 

isA/  =  l.  The  rotational  transitions  will  be  written  as  P{J lm,er ) and R(JLower ) .  The 

notation  used  for  vibrational  transitions  is  v(0,3)  which  means  the  transition  from  the 
ground  vibrational  state  to  the  3rd  vibrational  level  or  the  2nd  overtone,  and  v(3,2) 
describes  the  transition  from  the  3rd  vibrational  level  to  the  2nd  vibrational  level.  For 
example  consider  the  v(0,3)R(6)  transition,  this  rovibrational  transition  is  from  the  J=6 
rotational  level  within  the  ground  state  vibrational  manifold  to  the  3rd  vibrational 
manifold  at  the  J=7  rotational  level. 
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Table  2.1.  Spectroscopic  constants 


Constant 

Value  (cm"1) 

2169.8124 

<°eXe 

13.2879 

2. 

1.9312 

De 

6.12E-06 

ae 

0.0175 

Anharmonic  Potential  and  Other  Higher  Order  Terms 

The  above  results  for  a  diatomic  molecule  using  the  rigid  rotor  and  harmonic 
oscillator  models  provides  a  good  first  order  approximation,  but  a  more  precise  model  is 
required  to  adequately  model  the  second  overtone  pumping  of  CO.  Figure  2.3  shows  the 
harmonic  oscillator  is  good  approximation  for  small  ro vibrational  transitions,  Av  =  ±1 . 
However,  for  higher  overtone  transitions,  the  deviation  between  the  harmonic  and 
anharmonic  potential  grows  as  seen  in  Figure  2.3. 
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CO  Potential  Energy  Harmonic  &  Anharmonic 


Figure  2.3.  Harmonic  and  Anharmonic  potential  wells.  For  v  =  0,1  the  harmonic 
potential  is  a  close  approximation,  but  for  v  =  2+  it  begins  to  deviate. 

The  model  can  be  drastically  improved  by  including  the  anharmonic  oscillator  terms, 
rotational-vibrational  interactions,  and  rotational  effects.  The  harmonic  oscillator  model 
can  be  improved  by  including  higher  order  terms  in  the  Taylor  series  expansion  of  U(r) 
as  shown  in  equation  2.10. 

Rotational  transitions  are  super  imposed  upon  vibrational  transitions.  Therefore, 
it  is  convenient  to  write  the  equation  for  a  rovibrational  transition  [13]. 


Kj  -Ey+Ej  -  ae  (V  +  X)  +  Bj(j  +  l) 

(2.12) 

~  h  r  1 1 

(2.13) 

B  =  — - —  [cm  1 

8  7t2cI  L  J 
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Molecules  are  not  100%  rigid.  As  the  rotational  energy  grows  (higher  J  levels),  the 
centrifugal  force  causes  the  bonds  to  stretch  [14].  The  stretching  is  a  small  effect  and  can 
be  corrected  for  using  perturbation  theory  resulting  in  a  centrifugal  distortion  term,  D . 

Ej=BJ(J  +  \)-DJ2(J  +  \)2  (2.14) 

A  similar  effect  on  bond  length  occurs  with  vibrational  transitions,  causing  the  rotational 
energy  to  have  a  vibrational  dependency.  As  the  vibrational  energy  increases,  the  inter- 
nuclear  bond  length  also  increases  vibration-rotation  interaction  [14].  The  vibrational 
dependent  rotational  energy  is: 


Ej  =  BvJ(J  +  l)-DJ2(j  +  l )2 

(2.15) 

(2.16) 

The  harmonic  oscillator  model  must  also  be  improved  by  including  higher  order 
terms.  The  higher  order  terms  are  termed  anharmonicity  constants.  The  resulting 
potential  function  can  be  seen  in  Figure  2.2.  The  equation  for  the  vibrational  energy 
including  the  higher  order  terms  is: 

Ev  =6)e(v  +  y2)-coeZe{v  +  y2)2  +®eye(v+l/2y  (2.17) 

where  %e  « 1  •  Now  the  total  rovibronic  energy  of  a  diatomic  molecule  can  be  written 
as: 

SvJ=v„+ffl,(v+X)-ffl,Z.(v  +  K),+ffl.>'.(v  +  X)I+W  +  l)-Dy!p  +  0  (2-18) 

The  energy  associated  with  vibrational  transitions  is  approximately  100  to  1000  times 
that  of  rotational  transitions.  The  energies  and  Einstein  A  coefficients  associated  with 
common  transitions  in  this  research  project  can  be  found  in  Tables  2.2-2. 3. 
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Table  2.2.  The  energies,  wavelengths,  and  A  coefficients  for  v  (3,2). 


v  =  3  to  v  =2 

Jlower 

R-Branch 

P-Branch 

cmA-1 

nm 

A  (Hz) 

cmA-1 

nm 

A  (Hz) 

2094.1099 

4775.2986 

33.2477 

- 

- 

- 

i 

2097.8146 

4766.8655 

40.1147 

2086.5951 

4792.4966 

98.6617 

2 

2101.4840 

4758.5420 

43.2121 

2082.7854 

4801.2628 

65.4080 

3 

2105.1180 

4750.3275 

45.0525 

2078.9410 

4810.1414 

58.5366 

4 

2108.7164 

4742.2213 

46.3210 

2075.0620 

4819.1332 

55.4334 

5 

2112.2791 

4734.2229 

47.2810 

2071.1486 

4828.2388 

53.5859 

6 

2115.8059 

4726.3315 

48.0555 

2067.2009 

4837.4592 

52.3088 

7 

2119.2966 

4718.5467 

48.7097 

2063.2191 

4846.7949 

51.3386 

8 

2122.7512 

4710.8677 

49.2811 

2059.2034 

4856.2468 

50.5523 

9 

2126.1694 

4703.2941 

49.7930 

2055.1538 

4865.8158 

49.8848 

10 

2129.5512 

4695.8252 

50.2607 

2051.0706 

4875.5026 

49.2986 

11 

2132.8963 

4688.4605 

50.6944 

2046.9538 

4885.3080 

48.7702 

12 

2136.2047 

4681.1994 

51.1014 

2042.8037 

4895.2329 

48.2844 

13 

2139.4761 

4674.0415 

51.4870 

2038.6203 

4905.2782 

47.8310 

14 

2142.7105 

4666.9861 

51.8549 

2034.4039 

4915.4448 

47.4027 

15 

2145.9076 

4660.0328 

52.2082 

2030.1545 

4925.7334 

46.9943 

16 

2149.0674 

4653.1811 

52.5490 

2025.8724 

4936.1451 

46.6018 

17 

2152.1897 

4646.4305 

52.8792 

2021.5576 

4946.6808 

46.2224 

18 

2155.2744 

4639.7805 

53.2001 

2017.2103 

4957.3414 

45.8537 

19 

2158.3212 

4633.2307 

53.5128 

2012.8306 

4968.1279 

45.4941 

20 

2161.3301 

4626.7806 

53.8182 

2008.4188 

4979.0412 

45.1420 

21 

2164.3009 

4620.4297 

54.1171 

2003.9749 

4990.0824 

44.7965 

22 

2167.2334 

4614.1777 

54.4100 

1999.4991 

5001.2525 

44.4565 

23 

2170.1275 

4608.0241 

54.6974 

1994.9916 

5012.5524 

44.1214 

24 

2172.9831 

4601.9685 

54.9798 

1990.4525 

5023.9833 

43.7905 

25 

2175.8000 

4596.0106 

55.2574 

1985.8819 

5035.5462 

43.4633 

26 

2178.5781 

4590.1499 

55.5306 

1981.2800 

5047.2423 

43.1394 

27 

2181.3172 

4584.3860 

55.7997 

1976.6469 

5059.0725 

42.8185 

28 

2184.0171 

4578.7187 

56.0647 

1971.9828 

5071.0381 

42.5002 

29 

2186.6778 

4573.1474 

56.3259 

1967.2879 

5083.1402 

42.1843 

2.2  Selection  Rules 

The  rovibrational  selection  rules  ( Av  =  ±1 ,  A/  =  ±1)  stated  in  the  previous 
section  are  derived  from  the  time-dependent  version  of  the  Schrodinger’s  equation  [15]. 
The  time  dependent  Schrodinger’s  equations  is: 


dt 


(2.19) 


12 


(2.20) 


^(r,t)=^n(ryE^ 

where  equation  2.20  is  the  time  dependent  wave  equation,  and  ^(r)  refers  to  the 

stationary  states,  unchanging  with  time.  Now  consider  a  non-stationary  state  system 
where  an  electric  field  interacts  with  a  molecule.  The  electric  field  of  interest  will  be 
written  as: 

£  =  £  cos(2;rw)  (2.21) 

where  v  is  the  frequency.  The  Hamiltonian  operator  for  the  interaction  between  an 
electric  field  and  a  molecule  is  written  as: 

H®  =  -//•£  =  -/j.-E0cos{lmjt)  (2.22) 

where  the  superscript (1)  is  used  to  denote  the  perturbation  portion  of  the  Schrodinger’s 
equation  listed  below. 

H  =  H{0)+H{']  (2.23) 
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Table  2.3.  The  energies,  wavelengths,  and  A  coefficients  for  v  (2,1). 


v  =  2  to  v  =1 

Jlower 

R-Branch 

P-Branch 

cmA-1 

nm 

A  (Hz) 

cmA-1 

nm 

A  (Hz) 

hh 

2120.5647 

4715.7251 

23.0686 

- 

- 

- 

1 

2124.3044 

4707.4233 

27.8341 

2112.9799 

4732.6526 

68.4542 

2 

2128.0088 

4699.2286 

29.9844 

2109.1352 

4741.2798 

45.3804 

3 

2131.6778 

4691.1404 

31.2625 

2105.2558 

4750.0167 

40.6117 

4 

2135.3112 

4683.1581 

32.1440 

2101.3418 

4758.8641 

38.4576 

5 

2138.9089 

4675.2810 

32.8115 

2097.3934 

4767.8228 

37.1748 

6 

2142.4707 

4667.5085 

33.3503 

2093.4107 

4776.8935 

36.2878 

7 

2145.9964 

4659.8401 

33.8057 

2089.3939 

4786.0769 

35.6138 

8 

2149.4860 

4652.2751 

34.2037 

2085.3432 

4795.3737 

35.0675 

9 

2152.9392 

4644.8130 

34.5606 

2081.2586 

4804.7849 

34.6037 

10 

2156.3560 

4637.4533 

34.8867 

2077.1404 

4814.3111 

34.1962 

11 

2159.7361 

4630.1953 

35.1894 

2072.9886 

4823.9531 

33.8289 

12 

2163.0795 

4623.0387 

35.4736 

2068.8035 

4833.7118 

33.4913 

13 

2166.3859 

4615.9828 

35.7430 

2064.5851 

4843.5881 

33.1761 

14 

2169.6553 

4609.0271 

36.0002 

2060.3337 

4853.5827 

32.8784 

15 

2172.8874 

4602.1712 

36.2473 

2056.0493 

4863.6965 

32.5946 

16 

2176.0822 

4595.4146 

36.4858 

2051.7322 

4873.9305 

32.3218 

17 

2179.2395 

4588.7567 

36.7170 

2047.3824 

4884.2855 

32.0582 

18 

2182.3592 

4582.1972 

36.9418 

2043.0001 

4894.7624 

31.8020 

19 

2185.4410 

4575.7355 

37.1610 

2038.5854 

4905.3622 

31.5522 

20 

2188.4849 

4569.3713 

37.3752 

2034.1386 

4916.0859 

31.3076 

21 

2191.4907 

4563.1041 

37.5849 

2029.6597 

4926.9343 

31.0676 

22 

2194.4582 

4556.9335 

37.7906 

2025.1489 

4937.9084 

30.8316 

23 

2197.3873 

4550.8590 

37.9925 

2020.6064 

4949.0093 

30.5989 

24 

2200.2779 

4544.8804 

38.1910 

2016.0323 

4960.2381 

30.3692 

25 

2203.1298 

4538.9971 

38.3862 

2011.4267 

4971.5956 

30.1421 

26 

2205.9429 

4533.2089 

38.5785 

2006.7898 

4983.0830 

29.9173 

27 

2208.7170 

4527.5153 

38.7679 

2002.1217 

4994.7014 

29.6946 

28 

2211.4519 

4521.9161 

38.9546 

1997.4226 

5006.4518 

29.4738 

29 

2214.1476 

4516.4107 

39.1387 

1992.6927 

5018.3353 

29.2547 

The  H®  will  be  treated  as  a  small  perturbation  to  the  standard  time  dependent 


Hamiltonian,  equation  2.19.  The  selection  rules  will  be  derived  from  a  two-state  system 
in  which  the  wave  equation  is  written  as  a  linear  combination  of  two  wave  functions,  one 


for  each  state.  Each  state  is  multiplied  by  a  probability  function  that  depends  on  time.  At 
time  zero,  the  system  is  only  in  state  1 ,  a,  (t  -  0)  =  1 ,  a2  (t  =  0)  =  0  . 

Y(/)  =  a,(/>F1(/)+a2(^2(/)  (2.24)  ‘ 
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After  substituting  this  back  into  equation  2.23  and  after  some  simplification  and 
manipulation  the  following  results  are  found: 

°c  C  •  JvF2*//zvP1e?z'  (2.25) 

The  time  rate  of  change  of  ai  represents  the  rate  at  which  the  system  is  changing  from 
state  1  to  state  2.  C  is  simply  a  constant  term  containing  the  energies  of  the  two  states 
and  the  energy  of  the  incident  electric  field  and  does  not  determine  the  probability  of  the 
transition.  The  integral  in  equation  2.25  is  the  key  component  in  deriving  the  selection 
rules  and  represents  the  transition  dipole  moment.  It  is  obvious  the  transition  moment 
must  be  nonzero  in  order  for  a  transition  to  occur. 

Let  us  take  a  closer  look  at  the  transition  moment  function. 

Mu,  =  (2.26) 

This  equation  describes  how  “easy”  the  system  can  make  a  transition  from  state  1  to  state 
2,  and  is  the  key  to  determining  selection  rules  [15].  The  rotational  level  selection  rules 
will  be  derived  first.  Assume  the  electric  field  lies  only  along  the  z-axis,  and  simply 
substitute  equation  2.23  (/jz  =  //cos(<9))  and  the  solutions  of  the  rigid  rotor  model  into 
equation  2.26. 

2  nn 

(/L  )u  =  M  {  jyf  (0, </>)'  cos(0)  7/  (0, </>)sm(&)d6tf</>  (2.27) 

0  0 

where  the  dipole  moment  is  pulled  outside  the  integral,  which  directly  shows  there  must 
be  a  permanent  dipole  moment  present  in  order  for  the  molecule  to  have  a  pure  rotational 
spectrum  [15].  The  solution  to  the  integral  in  equation  2.27  is  only  non-zero  for  values 
in  which  J'  =  J  ±1 . 
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The  vibrational  transition  selection  rule  follows  a  similar  path.  It  was  stated 
above  that  only  Av  =  ±l  transitions  are  allowed.  This  is  true  for  a  first  order 
approximation.  First,  recall  the  solutions  to  the  nuclear  portion  of  the  harmonic-oscillator 
model  are  Hermite  polynomials  [15]. 

VM)  =  NvHv{ay>qyaq^  (2.28) 

After  applying  the  same  substitutions  as  before  in  equation  2.26  the  following  results: 

co 

=  \n,h, 

—00 

In  order  to  solve  the  integral,  the  dipole  moment  must  be  expanded  around  the 
equilibrium  nuclear  separation.  Including  only  the  first-order  term  results  in  the  Av  =  ±1 
selection  rule.  However,  when  the  higher  order  terms  are  included  in  the  dipole  moment 
expansion,  the  allowed  vibrational  transitions  are  expanded  to  Av  =  ±1,±2,±3...  depending 
on  how  many  nonlinear  terms  of  the  dipole  moment  expansion  are  included.  The  higher 
order  terms  are  much  smaller  than  lower  order  terms.  This  means  that  the  bigger  the  Av 
transition  the  more  unlikely  it  is  to  occur  [13].  In  CO  a  Av  =  ±2  is  approximately  150 
times  more  likely  than  a  Av  =  ±3  .  All  Av  >  1  re  called  overtone  transitions. 

2,3  Thermal  Population 

In  the  previous  theory  section,  the  translational  energy  of  the  diatomic  molecule 
was  excluded.  It  was  shown  that  the  energy  of  a  diatomic  molecule  can  be  written  as  a 
sum  of  the  individual  energies:  electronic,  vibrational,  rotational,  and  now  translational 

[14]. 


a^qy  '2  '2  dq  (2.29) 
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£Tot  =  £ elec  +  £vi  b  +  £rol  +  £  Irons  (2.30) 

The  translational  energy  is  of  particular  interest  because  physical  collisions  between 
molecules  cause  quantum  mechanical  transitions.  In  other  words,  the  collision  between 
two  diatomic  molecules  can  cause  a  shift  in  the  molecules’  rotational  level,  inducing  a 
rotational  transition.  These  collisions  which  are  responsible  for  the  molecules’  thermal, 
distribution  are  described  by  the  Boltzmannn  relationship, 


F.  =■ 


Sse 


zF 

kT 


Hsj' 


-II 

kT 


(2.31) 


where  Fs  is  the  fraction  of  molecules,  Es  is  the  energy,  gs  is  the  degeneracy  of  state  s,  and 
the  denominator  is  the  summation  over  all  energies  [13].  Notice  the  Boltzmannn 
relationship  is  a  function  of  the  Boltzmannn  constant  times  the  temperature,  which  is 
approximately  190  cm1  at  room  temperature.  The  approximate  energy  of  a  diatomic 
molecule  associated  with  a  vibrational  transition  is  on  the  order  of  2000  cm'1.  This 
shows  that  almost  all  of  the  molecules  at  room  temperature  will  be  in  the  ground  state. 
Rotational  energy  levels  are  spaced  much  more  finely  than  vibrational  levels,  implying  at 
room  temperature  there  will  be  a  thermal  distribution  of  molecules  throughout  many 
rotational  levels,  Figure2.4. 
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Carbon  Monoxide 
Thermal  Vibrational  Distribution 


Vibrational  Level 
Carbon  Monoxide 


Figure  2.4.  Vibrational  and  rotational  distribution  of  CO  at  room  temperature. 

2.4  Kinetic  Processes 

This  section  will  cover  the  various  kinetic  processes  associated  with  optically 
pumped  CO  lasers.  Rotational  relaxation  is  the  only  process  that  occurs  on  the  same  time 
scales  as  lasing.  Vibrational-vibrational  and  vibrational-translational  transitions  occur  on 
time  scales  much  larger  than  the  lasing  processes.  It  will  be  shown  in  the  numerical 
model  chapter  that  V-V  and  V-T  processes  can  be  neglected  for  these  experiments. 
Rotational  relaxation  is  key  to  the  scaling  and  efficiency  of  the  CO  optically  pumped 
laser.  Its  major  role  is  discussed  in  the  radiative  process  section. 
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Vibrational-Vibrational  (V-V)  Transitions 

V-V  transitions  occur  when  two  molecules  transfer  vibrational  energy  among  one 
another.  Typically,  one  molecule  transitions  up  a  vibrational  level  and  the  other  drops 
down  a  vibrational  level. 

4v)  +  5(v')-j£->4v-l)  +  5(v'  +  l)  +  A£:  (2.32) 

Using  the  calculations  from  McCord  (1999),  the  rate  constant  for  V-V  relaxation  at  room 
temperature  is  kf  =  2.73  •  1(T12  cmYmo,e.s .  This  leads  to  a  V-V  relaxation  of  630  ns  with  20 

Torr  of  CO  [4],  This  is  well  beyond  the  lasing  time  scales  of  the  CO  cascade  laser,  which 
are  on  the  order  of  50-100  ns.  The  effect  of  V-V  transitions  in  the  numerical  rate 
equation  model  will  be  shown  to  be  negligible. 

Vibrational-Translational  (V-T)  Transitions 

V-T  transitions  imply  a  molecule  imparts  its  vibrational  energy  to  a  collision 
partner  in  the  form  of  translational  energy.  V-T  is  an  extremely  slow  process,  and  with  a 
predicted  relaxation  time  of  60  ms  at  20  Torr  of  CO  [4],  Thus,  V-T  is  two  orders  of 
magnitude  slower  than  the  radiative  processes  in  CO  lasing  and  can  be  neglected. 

Rotational-Relaxation  (R-R) 

R-R  relaxations  occur  within  the  same  vibrational  manifold  when  two  molecules 
collide.  This  results  in  a  very  small  energy  release  and  can  be  written  as: 

CO(v,  J)  +  M  — kj—>  CO(v,  J')  +  M  +  AE  (2.33) 
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R-R  relaxations  occur  on  the  same  time  scale  as  the  CO  radiative  processes  and  must  be 
accounted  for  in  the  rate  equation  model.  The  time  scales  of  R-R  will  be  explored  in 
greater  detail  in  a  following  section. 

A  qualitative  discussion  of  R-R  effects  on  the  lasing  properties  of  an  optically 
pumped  CO  laser  will  follow.  First,  let  us  consider  two  cases  of  the  R-R  rates.  Case  1 
will  represent  near  instantaneous  R-R,  and  case  2  will  assume  a  relatively  slow  R-R  rate. 

Absorption  with  Instantaneous  R-R  and  Slow  R-R  Conditions 

The  near  instantaneous,  1  ns,  R-R  case  implies  that  a  vibrational  manifold  will 
instantaneously  reach  the  thermal  Boltzmannn  distribution.  The  first  benefit  of 
instantaneous  R-R  is  higher  pump  pulse  absorption  [17].  This  is  obvious  by  simply 
examining  the  laser  gain/absorption  equation. 

r.„  («)  =  o* (uf  JV„  -  ■ iiL N' )  (2.34) 

{  Sl  J 

This  equation  will  be  explained  in  much  greater  detail  in  the  radiative  process  section, 
but  for  now  we  will  quantitatively  examine  it.  The  crabs  (o)  is  the  absorption  cross 
section,  and  can  be  treated  as  a  constant  for  now.  N(J ,  NL  are  the  populations  of  the 

upper  and  lower  states  respectively,  and  —  is  simply  a  degeneracy  term  and  can  be 

Sl 

assumed  to  be  one  for  this  discussion.  The  upper  and  lower  states  in  this  equation  are 
rovibrational  levels;  they  are  specific  rotational  levels  within  a  specific  vibrational 
manifold.  yahs{°)  represents  the  absorption  per  centimeter  of  the  incoming  pump  pulse. 
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Notice  once  the  upper  state,  Nv ,  reaches  lA  the  population  of  the  lower  state,  NL ,  the 

pump  pulse  absorption  goes  to  zero.  This  is  termed  saturation. 

Now  consider  the  case  of  slow  R-R.  This  implies  the  pump  pulse  will  remove 
molecules  from  a  specific  rotational  level  in  the  lower  state  and  place  them  in  a  specific 
rotational  level  in  a  higher  vibrational  state.  The  population  inversion  will  remain 
constant,  and  the  transition  will  reach  saturation  quickly,  since  only  the  molecules  in  a 
specific  rovibrational  level  can  absorb  the  pump  pulse. 

Now  imagine  if  the  R-R  process  happened  instantaneously.  The  lower  rotational 
level  would  fill  in  from  adjacent  rotational  levels  as  quickly  as  the  pump  removed  them. 
This  means  NL  can  now  be  considered  the  population  of  the  entire  lower  vibrational 
manifold.  The  upper  rotational  level  will  relax  out  into  adjacent  states,  spreading  its 
population  throughout  the  entire  vibrational  manifold  in  effect  lowering  Nv .  This 

means  the  pump  pulse  can  now  pull  CO  molecules  from  the  entire  lower  vibrational 
manifold  instead  of  a  specific  rovibrational  level,  and  the  upper  rovibrational  level  never 
fills  up.  This  allows  for  higher  pump  absorption,  and  greater  efficiencies.  One  may 
question  that  if  R-R  is  instantaneous  and  if  the  population  inversion  spreads  out  across 
the  entire  vibrational  level,  the  gain  must  also  decrease.  This  is  true;  however,  the 
decreased  gain  is  actually  another  benefit  of  instantaneous  R-R. 

Determining  R-R  Rate  Constants 

Phipps  made  the  first  direct  measurements  of  CO-CO  rotational  relaxation 
constants  using  a  double  IR-IR  pump  probe  resonance  technique.  The  results  of  Phipps’s 
research  showed  the  statistical  exponential  power  gap  (SPEG)  law  produced  the  most 
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accurate  rate  constants  for  CO-CO  interactions  [18].  The  master  equation  governing  the 
evolution  of  the  rotational  population  in  a  specific  v  level  is  [19] 

(2.35) 

at  f 

where  Nj  is  the  population  of  the  Ith  level,  and  kJ(_f  is  the  rate  constant  for  the  transition 

from  the  J’  to  the  J  level  [19].  Assuming  the  thermal  distribution  of  CO  at  room 
temperature  only  populates  the  first  twenty  rotational  levels  equation,  2.35  would  require 
400  rate  constants  to  completely  describe  the  system,  and  no  experimental  results  contain 
enough  information  to  determine  each  rate  constant  independently.  The  approach  taken 
by  Phipps  is  to  solve  the  coupled  differential  equations  2.35  and  compare  the  solution  to 
the  measured  data.  The  rate  constants  used  in  the  numerical  solution  to  equation  2.35  are 
generated  from  various  state  to  state  fitting  laws.  The  SPEG  law  produced  the  best 
results  with  the  least  amount  of  error.  The  SPEG  law  is  defined  as 


,  [  A Erj 

kf<-j  -  ai 


(2.36) 


where  Z?v=3  is  the  CO  rotational  constant,  A EfJ  is  the  energy  gap  of  the  transition,  and  the 

other  coefficients  are  fitting  parameters.  The  best  fit  parameters  for  the  SPEG  law 
determined  by  Phipps  [19]  where  used  in  the  rate  equation  model  described  later. 
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2.5  Radiative  Processes 


The  three  radiative  processes  developed  by  Einstein  include  spontaneous 
emission,  absorption,  and  stimulated  emission.  All  three  processes  describe  a 
relationship  between  light  and  the  population  of  the  system’s  state  [20].  Each  process 
will  be  discussed  in  detail. 

Spontaneous  emission  is  the  process  by  which  an  atom  from  state  2  spontaneously 
decays  to  state  1  with  in  a  vacuum  field.  This,  in  turn,  states  that  a  system  with  an 
excited  population  will  decay  to  the  ground  state  through  only  spontaneous  emission  at  a 
rate  of  A21~' .  The  spontaneous  emission  can  be  written  as  a  simple  rate  equation. 


dN, 


dt 


Spon  tan  eous 
Emission 


~  ^21-^2 


(2.38) 


The  energy  must  be  conserved;  therefore  when  the  atom  decays,  a  photon  is  released  with 
the  corresponding  energy.  The  released  photon  is  of  random  direction,  phase,  and 
polarization. 

Absorption  is  the  process  in  which  an  atom  from  state  1  absorbs  a  photon  and 
excites  the  atom  to  the  upper  state  2.  The  photon  energy  must  match  the  atomic 
transition  energy: 


dN2 

dt 


Absorption 


BaN,p(y)  =  -^l 


dt 


Absorption 


(2.39) 


where  p{o)  is  the  energy  density  per  frequency  interval.  Notice  the  relation  between  the 
upper  and  lower  states.  As  the  lower  state  looses  an  atom  the  upper  state  gains  an  atom; 
and  total  population  remains  constant. 
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Stimulated  emission  is  the  main  process  within  a  laser  medium.  Stimulated 
emission  is  the  result  of  a  photon,  generated  from  a  transition,  colliding  with  an  atom  in 
state  2,  and  causing  the  atom  to  drop  to  state  1  releasing  another  photon. 


dN, 


dt 


„  ,T  .  .  dN, 
=  -B2  lNlp(v)  =  — -j- 

Stimluated  Qt 

Emission 


Stimulated 

Emission 


(2.40) 


The  photon  generated  from  stimulated  emission  must  have  the  exact  same  frequency, 
polarization,  direction,  and  phase  of  the  incoming  photon  [20]. 

Einstein  combined  the  three  radiative  processes  under  the  framework  of 
thermodynamic  equilibrium  to  develop  key  relationships  between  the  coefficients. 
Basically,  the  combination  of  the  three  rate  equations  describes  the  system’s  population 
state  at  any  given  time. 


dN2 

dt 


Radiative 


-A2lN2  +  BX2Nxp(v)-B2XN2p(v)  =  - 


dNx 

dt 


Radiative 


(2.41) 


This  equation  shows  as  state  2  loses  an  atom,  state  1  gains  an  atom.  Assuming 
equilibrium  conditions  and  applying  Boltzmann  statistics,  Einstein  showed  the  following 
relationships  between  the  radiative  coefficients 

§2^2]  ~ 

A2 ,  _  8 nhv* 

B2]  c 

where  gx  2  is  the  degeneracy  of  states  1  and  2  respectively.  The  above  relationships  show 
that  if  one  process  is  known,  the  others  are  also  known  [21]. 


(2.42) 

(2.43) 
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2.6  Line  Shape,  Gain,  and  Absorption 

The  radiative  processes  involve  photons  of  specific  energies: 

hv  =  E2-Ex  (2.44) 

where  E2 ,  are  the  energies  of  the  states  2  and  1  respectively.  This  implies  an  exact 

knowledge  of  the  energy  of  an  atom,  which  is  impossible  due  to  the  uncertainty  principle. 
This  uncertainty  results  in  a  bandwidth  or  a  distribution  of  energies  around  a  central 
value  when  many  atoms  are  considered  in  the  system.  This  distribution  is  called  the  line 
shape,  and  its  shape  is  influenced  by  many  parameters. 

Line  Shape 

The  line  shape  function  is  defined  as  g(v),  and  g(y)dv  is  the  probability  a  photon 
from  the  radiative  processes  will  appear  in  the  interval  vto  v  +  dv  [21],  The  same  line 
shape  function  applies  for  all  radiative  processes.  The  line  shape  is  considered  broadened 
when  the  spectral  width,  Av,  increases.  There  are  two  types  of  broadening, 
homogeneous  and  inhomogeneous. 

Pressure  Broadening 

Homogeneous  broadening  is  caused  by  a  mechanism  that  affects  all  atoms  equally  [21]. 
Pressure  broadening  is  the  most  important  homogenous  broadening  effect  and  must  be 
accounted  for.  Pressure  broadening  is  a  result  of  two  molecules  elastically  colliding, 
which  results  in  a  discontinuous  phase  jump  in  the  radiation  [21].  The  phase  jump 
appears  as  broadening  in  frequency  space.  The  collision  rate  of  the  molecules  directly 
influences  the  pressure  broadening  effect.  The  more  CO  molecules  that  are  present 
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(higher  pressure),  implies  more  collisions  which  in  turn  broadens  the  line  more.  The 
same  argument  can  be  made  for  an  increase  in  temperature  and  the  addition  of  different 
gas  species  or  collision  partners.  The  FWHM  of  the  pressure  broadening  is  given  by: 


5vP=Y,ai(r)Pi 


(2.45) 


where  cr,(r)  is  the  temperature  dependent  pressure  broadening  coefficient  for  species  i 
[MHz/Torr]  and  p,  is  the  partial  pressure  of  species  i.  Pressure  broadening,  collision 
broadening,  is  represented  by  a  Lorentzian  line  shape  function, 


/(v)  = 


Awl 


In 


(  \2  WL 

(v'-O  +  ~r 


(2.46) 


wL  = 


NQ  1 6kT 
n  V  n. M 


(2.47) 


where  wL  is  the  Lorentzian  width,  vc  is  the  center  frequency,  A  is  the  spectral  area,  N  is 


the  number  density  of  the  state,  M  is  the  reduced  mass  and  Q  is  the  collisional  cross 
section  [13]. 


Table  2.4.  CO  self-pressure  broadening  coefficients  [22]. 


R(J) 

Wavelength 

(nm) 

Broadening 

(MHz/mBar) 

0 

1573.7693 

2.556 

i 

1572.8695 

2.328 

2 

1571.9969 

2.244 

3 

1571.1512 

2.213 

4 

1570.3325 

2.099 

5 

1569.5406 

2.04 

6 

1568.7757 

1.985 

7 

1568.0375 

1.948 

8 

1567.3261 

1.912 

9 

1566.6415 

1.871 

10 

1565.9836 

1.854 

R(J) 

Wavelength 

(nm) 

Broadening 

(MHz/mBar) 

11 

1565.3523 

1.835 

12 

1564.7477 

1.802 

13 

1564.1698 

1.767 

14 

1563.6184 

1.742 

15 

1563.0936 

1.715 

16 

1562.5954 

1.664 

17 

1562.1237 

1.682 

18 

1561.6786 

1.637 

19 

1561.2600 

1.591 

20 

1560.8679 

1.571 
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Doppler  Broadening 

Inhomogeneous  broadening  pertains  to  the  line  broadening  of  different  groups  of 
atoms  that  are  broadened  differently.  Doppler  broadening  is  an  inhomogeneously 
broadening  process  since  it  depends  on  the  velocity  of  the  gas  flow.  The  velocity 
distribution  can  be  represented  by  a  Maxwell  Boltzmann  distribution,  which  in  turn 
yields  a  Gaussian  line  shape, 


1 41n(2) 

exp 

f - 

o 

1 

> 

CN 

1 

1 _ 

2" 

V  n 

(2.48) 


where  wG  is  the  Gaussian  line  width,  vc  is  the  center  frequency,  M  is  the  molecular  mass, 
T  is  the  gas  temperature,  and  R  is  the  ideal  gas  constant  [13]. 


Voigt  Profile 

At  low  pressures,  1-5  Torr,  the  line  shape  of  CO  can  be  represented  by  using  a 
Gaussian  profile.  However,  the  majority  of  the  laser  experiments  are  at  higher  CO 
pressures,  therefore  pressure  broadening  can  not  be  neglected.  This  forces  the  line  shape 
to  be  a  convolution  of  the  Gaussian  and  Lorentzian  line  shape:  a  Voigt  function.  The 
Voigt  function  is: 


,(  \  .  21n2  wL  °r 


e"  dt 


(Vtaim  +[(/4h2) 


V-  v. 


\2 


(2.50) 


wr. 


All  the  variables  are  the  same  as  the  Gaussian  and  Lorentzian  line  shapes.  All  line 
shapes  mentioned  above  can  be  seen  in  Figure  2.5. 
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Figure  2.5.  Lorentzian,  Gaussian,  and  Voigt  profiles.  The  Lorentzian  profile  (green) 
represents  only  pressure  broadening,  the  Gaussian  profile  (blue)  show  only  Doppler 
broadening  effects,  and  the  Voigt  profile  (red)  is  a  convolution  of  both  pressure 
broadening  and  Doppler  broadening. 


Gain 


Assuming  the  radiant  energy  density  spectral  width  is  much  smaller  than  the 
spectral  width  of  g(v) ,  the  radiative  rate  equation  2.41  can  be  written  as: 


dN, 


dt 


—  A2XN2 


°{v)h 


Radiative 


ho 


N2  ~  —  Nx 

8i 


(2.51) 


°{v)=a2X  rsM 

8m 


(2.52) 


where  — ^  is  the  photon  flux  of  the  lasing  transition  and  cr(v)  is  the  stimulated  emission 
hv 

cross  section  with  units  of  area.  Typically  the  first  term  in  equation  2.51,  spontaneous 
emission,  is  dropped  since  it  is  very  small  compared  to  the  stimulated  emission. 
However,  the  spontaneous  emission  is  very  important  since  it  is  the  starting  process  of  a 
laser.  The  cross  sectional  area  is  the  area  of  the  atom  visible  to  the  photon  flux  [21]. 
Notice  the  absorption  cross  section  depends  directly  on  the  Einstein  A-coefficient.  The 
A-coefficient  is  derived  from  a  quantum  mechanical  approach,  and  it  is  directly  related  to 
the  transition  dipole  moment  which  was  discussed  in  section  2.2.  As  the  transition 
becomes  less  likely  (more  forbidden)  the  transition  dipole  becomes  smaller,  resulting  in  a 
smaller  cross  sectional  area.  This  simple  statement  defines  the  largest  obstacle  in  the  2nd 
overtone  optically  pumped  CO  laser.  Since  the  pumped  transition  involves  Av  =  2 ,  it 
has  a  very  small  cross  section.  The  small  radiation  absorption  by  the  CO  at  the  Av  =  2 
wavelength  is  a  direct  result  of  the  small  cross  section.  Another  interesting  result  of 

ar 

equation  2.51  is  the  fact  that  once  the  population  of  the  upper  state  is  equal  to  —Nx,  the 

g\ 

stimulated  emission  stops,  in  other  words  the  laser’s  gain  saturates. 
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Neglecting  the  spontaneous  emission  term  and  multiplying  by  the  photon  energy, 
ho  equation  2.44  can  be  written  as: 


*v"> 


dt 


=  -°(v) 


Radiative 


N2-^-Nl 
g  l 


(2.53) 


The  term  on  the  left  side  represents  the  energy  of  a  released  photon  times  the  rate  of 
decay,  and  the  right  side  term  represents  the  changing  intensity  of  the  emitted  photons. 
Since  the  photons  released  via  stimulated  emission  are  identical,  the  right  hand  side  of 
equation  2.53  is  the  change  of  the  radiated  intensity  with  respect  to  distance  through  the 
system.  The  following  results, 


—  JV, 
g\ 


=  r{v)iu 


(2.54) 


where  y{v)  is  the  gain  coefficient  of  the  system,  and  I  is  the  lasing  intensity.  It  is 
obvious  from  equation  2.54  that  a  population  inversion  must  be  created  in  order  for  the 
optical  system  to  have  gain.  The  gain  coefficient  is  written  as: 


r(v)  =  cr(v 


N2-^-Nl 

g\ 


(2.55) 


Absorption:  Beer-Lambert 

Solving  the  differential  equation  2.54  by  integrating,  results  in  the  well  known 
Beer-Lambert  law. 

(2.56) 

where  Io  is  the  initial  intensity  and  I  is  the  transmitted  intensity.  Replacing  the  gain 
coefficient  in  equation  2.56  results  in  the  following: 


y{v)  =  j\n\ 
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(2.57) 


f  g2l 

CT(V)  in 

"7°1 

UJ 

L 

I 

N,= 


This  simple  equation  states  the  number  density  of  the  lower  state  can  be  determined  by 
measuring  the  absorption  of  incident  radiation  through  a  medium,  where  I0  is  the  initial 
radiation  intensity  (before  the  medium)  and  I  is  the  output  intensity. 


Amplified  Spontaneous  Emission  (ASE) 

Amplified  spontaneous  emission  occurs  in  high  gain  laser  systems  as  a  result  of 
the  spontaneous  emission  being  amplified  by  the  same  group  of  atoms  [23].  The 
resulting  pulse  from  ASE  is  semi-coherent  and  contains  a  large  amount  of  energy.  An 
ASE  pulse  will  be  spectrally  coherent  and  have  a  moderate  divergence  which  is 
acceptable  for  some  applications.  However,  for  the  CO  cascade  laser,  a  highly  coherent 
output  beam  is  desired.  Thus,  ASE  poses  a  problem.  ASE  radiation  can  become  strong 
enough  to  cause  saturation  and  to  remove  a  significant  portion  of  the  energy  from  the 
resonator  by  pulling  atoms  from  the  excited  state  down,  which  also  results  in  a  lifetime 
shortening  [23].  The  maximum  gain  allowed  before  ASE  becomes  a  dominating  factor  is 
given  by  Verdeyen  as, 


G0(v0)  \  =  AnL 

A 

(2.58) 

(2.59) 

where  L  is  the  gain  length,  A  is  the  aperture  area,  and  G0  (i/0 )  is  the  power  gain  [24]. 
Equation  2.58  simplifies  into  an  expression  relating  the  gain  to  the  ratio  of  the  gain 
medium  radius  and  gain  length, 
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(2.60) 


.  (2L\ 

y0L  =  21n  — 

V  a  ) 

where  a  is  the  gain  diameter.  It  is  easily  seen  from  equations  2.58  and  2.59  that  ASE  is 
directly  dependent  on  the  gain  of  the  system  and  the  length  to  aperture  ratio. 

Laser  Gain  with  Instantaneous  R-R  and  Slow  R-R 

First,  assume  R-R  is  a  slow  process  compared  to  the  lasing  time  scales.  This 
creates  a  huge  population  inversion  between  the  pumped  rovibrational  level  and  the 
lower  lasing  level.  Since  all  the  inverted  population  goes  to  a  single  upper  rovibrational 
level  and  the  corresponding  lower  rovibrational  population  is  zero.  This  can  easily 
produce  gains  on  the  order  of  1000%  per  centimeter  [18],  Gain  at  such  extreme  levels 
produces  a  huge  ASE  pulse,  which  is  not  desired. 

This  is  where  instantaneous  R-R  comes  into  play.  If  the  R-R  is  instantaneous 
then  the  excited  rovibrational  level  will  spread  out  to  all  rotational  levels  within  the 
vibrational  manifold  based  on  the  Boltzmann  thermal  distribution.  This  in  turn 
distributes  the  gain  to  many  rovibrational  lines  throughout  the  entire  vibrational  manifold 
in  effect  lowering  the  gain  on  any  one  specific  lasing  line  to  more  reasonable  levels.  The 
laser  will  lase  from  the  entire  vibrational  manifold  from  many  rotational  levels  or  from  a 
single  rotational  level  depending  upon  the  resonator  setup.  All  without  ASE  losses. 

Through  a  quantitative  discussion  on  R-R  effects,  it  is  obvious  a  fast  R-R  rate  on 
order  of  ns  or  faster  is  desired.  This  provides  many  benefits  to  CO  optically  pumped 
lasers. 
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2.7  The  Resonator  and  Cavity  Modes 

A  simple  stable  curved-flat  resonator  was  used  for  all  experiments,  and  the 
theoretical  discussion  will  center  on  the  parameters  of  this  resonator.  The  laser  resonator 
provides  positive  feedback  to  the  amplification  loop,  allowing  a  laser  to  oscillate. 


Axial  Modes  of  a  Resonator 

Axial  modes  occur  in  a  passive  resonator.  They  result  from  the  phase  matching 
conditions  of  a  standing  wave  resonator  after  one  round  trip.  This  simply  means  the 
round  trip  length  must  be  an  integer  number  of  wavelengths  at  the  oscillating  frequency 
[25].  The  axial  mode  frequencies  are  easily  calculated  using  the  following  equation: 


vq=q 


2  L 


(2.61) 


where  q  is  an  integer,  and  vq  represents  a  set  of  frequencies  oscillating  in  the  optical 
cavity.  The  axial  mode  spacing  is  the  free  spectral  range  of  the  cavity. 

A  v=FSR  =  —  (2.62) 

q  2  L 

The  number  of  axial  modes  present  in  a  resonator  is  a  function  of  the  resonator’s  free 
spectral  range. 
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Transverse  Modes  of  a  Resonator 

Transverse  modes  are  the  radiation  intensity  distribution  in  a  cross-sectional  plane 
perpendicular  to  the  laser  axis  [26].  The  simplest  or  lowest  order  transverse  electric  and 
magnetic  mode  (TEM)  is  the  TEMoo.  The  TEMoo,  commonly  referred  to  as  the 
fundamental  mode,  and  has  the  intensity  profile  of  a  single  Gaussian  shaped  dot.  A  laser 
will  operate  in  all  modes  that  have  gain  above  the  threshold  level,  so  a  typical  laser  will 
oscillate  in  multiple  transverse  modes.  See  Figure  2.6  for  an  intensity  profile  of  higher 
order  modes. 


TEMoo  TEMoi  TEMn 


Figure  2.6.  Transverse  Mode  Patterns. 

Mode  Spacing 

A  laser  cavity  oscillates  on  all  the  frequencies  of  the  axial  and  transverse  modes 
with  gains  above  the  threshold  level.  This  typically  results  in  many  modes  operating 
simultaneously  and  independently  in  the  case  of  inhomogeneously  broadened  gain 
medium. 
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2.8  Mode  Locking 

Mode  locking  a  laser  is  defined  as  locking  the  phase  of  all  the  axial  modes 
together  at  a  specific  phase.  For  simplicity  assume  <f>n  =  0 .  A  very  simple  model  of  mode 
locking  will  be  laid  out  in  the  following  section;  the  model  will  predict  the  mode  locked 
pulse  width,  pulse  spacing,  and  peak  power. 

A  very  general  expression  for  the  total  electric  field  in  a  resonator  as  a  function  of 

time  is, 

£(0=Z£«cos(2;n/"/+^)  (2-65) 

n= 0 

where  En,  vn ,  and  (f>n  represent  the  amplitude,  frequency,  and  phase  of  the  nth  mode  at 
time  zero  [28],  E(t)  contains  the  individual  contributions  of  all  the  axial  modes. 
Equation  2.65  assumes  there  is  only  one  TEM0o  mode,  which  implies  vn  =v0+nAvq. 

Assuming  each  mode  has  the  same  amplitude  allows  En  to  be  removed  from  the 
summation.  This  is  obviously  a  gross  approximation  and  assumes  a  a  square  gain  profile 
versus  a  Voigt  profile.  The  time  dependent  intensity  can  be  found  by  squaring  equation 
2.58. 

i(‘)=h 

The  plot  of  the  intensity  versus  time  is  shown  in  Figure  2.8.  By  setting  the  phase  term 
equal  to  a  constant,  (j>n  =  0 ,  forces  a  mode-locked  condition  in  which  the  phase  of  each 
mode  is  the  same. 


AM 


n= 0 


i{2nnhvqt+<l>n) 


(2.66) 
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Time 


Figure  2.8.  Plot  of  4  modes  oscillating  with  random  phases. 


The  number  of  modes  allowed  to  oscillate  in  a  cavity  depends  on  a  few  factors.  First,  the 
resonator  length  determines  the  axial  mode  space  by  equation  2.53.  Each  mode  must  lie 
within  the  gain  profile,  such  that  the  maximum  number  of  oscillating  modes  is 

N  =  AFffTYAf//  ^  ^  where  A FWHM  is  the  full  width  half  maximum  of  the  gain  profile. 

This  assumes  all  the  modes  contained  within  the  gain  profile  are  above  threshold.  A 
square  gain  profile  is  not  a  good  approximation  of  the  gain  profile;  a  better 
approximation  is  a  Gaussian  profile.  In  Figure  2.7  it  is  easily  seen  not  all  the  modes 
under  the  Gaussian  profile  break  threshold;  thus  only  the  modes  within  the  FWHM  of  the 
line  shape  will  lase.  A  simple  expression  for  the  number  of  modes  assuming  a  Gaussian 


line  shape  is, 


A FWHM  ( ln(A~)y 

Av?  lln(2). 


(2.67) 
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where  K  is  simply  the  ratio  of  the  threshold  gain  and  the  small  signal  gain,  K  ~ 


Using  the  new  definition  of  the  number  of  modes  in  equation  2.67  and  assuming  a  locked 
phase,  the  following  results. 


I{t)  =  h 


sin2  (vV;rAvy) 
sin2(^Av?t) 


(2.68) 


As  more  modes  are  locked  together  the  more  intense  and  narrower  each  pulse  becomes. 
Notice  from  equation  2.68  there  are  two  ways  to  increase  the  number  of  modes  in  the 
laser  system.  First,  the  gain  profile,  A FWHM ,  can  be  broadened  allowing  more  modes  to 
fit  under  the  curve.  The  resonator  can  also  be  lengthened  to  decrease  A  vq .  The  FWHM 


of  the  mode-locked  pulse  is  given  by: 


A FWHM^%[2) 

Notice  by  increasing  tsFWHM ,  the  pulse  width  decreases;  however  by  decreasing 
Ar^the  pulse  width  does  not  change.  Figure  2.9  shows  a  pulse  train  of  4  and  8  mode 

locked  pulses. 
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Broadened  Gain 
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Figure  2.10  provides  a  rapid  visual  comparison  between  the  Fourier  transform  of  the  4 
and  8  mode  locked  pulse  trains.  The  beats  can  be  described  as  axial  mode  beats  and  are 
separated  by  A  uq. 

The  mode-locked  pulse  width  in  frequency  space  is  also  related  to  the  Finesse  of 
the  cavity  by: 


F  = 


»(*,*J* 

i -(*,*,)* 


(2.70) 


where  F  is  the  cavity  Finesse  and  the  R’s  are  the  reflectivities  of  the  resonator  mirrors. 


A  oT  [29], 


CHAPTER  THREE:  EXPERIMENTAL  AND  PROCEDURES 


This  chapter  outlines  in  detail  all  experimental  components  and  procedures  used 
to  collect  accurate  and  precise  CO  laser  pulse  data.  The  experimental  setup  is  shown 
schematically  in  Figure  3.1.  The  OPO  provided  a  tunable  single  mode  pump  pulse 
wavelength  matched  to  specific  rovibrational  transitions  in  the  2nd  overtone  of  CO. 
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Figure  3.1.  The  top  level  experimental  setup.  The  red  line  represents  the  1.5  micron 
pump  beam.  The  green  line  is  the  4.7  micron  CO  laser  pulse.  The  above  setup  provides 
a  double  pass  pump  run. 
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3.1  OPO  Specifics 


The  optical  parametric  oscillator  (OPO)  was  designed  and  manufactured  by 
Aculight  Corporation  through  a  contract  with  the  USAF.  The  OPO  was  used  as  the 
laser’s  pump  source,  providing  8ns,  10-20mJ  pulses  tunable  between  1.55  and  1.58 
microns.  The  OPO  was  wavelength  tuned  to  match  a  specific  rovibrational  transition  in 
the  CO  gas.  A  detailed  schematic  of  the  OPO  is  found  in  Figure  3.2. 


CCD  Camero 


Figure  3.2.  The  OPO’s  setup.  The  black  line  is  the  1064  nm  seeded  Nd:YAQ  pulse. 

The  red  line  is  the  single  mode  tunable  (1550-1580nm)  continuous  seed  laser  beam.  The 
orange  line  is  the  narrow  bandwidth  pump  pulse  output  from  the  OPO. 


42 


Nd:YAG 


A  Continuum  Q-switched  Nd:YAG  provided  the  energy  for  the  OPO  resonator. 
The  Nd:YAG  was  an  injection  seeded  YAG,  which  allowed  it  to  lase  on  a  single  axial 
mode.  See  Figure  3.3  for  the  seeded  versus  unseeded  Nd:YAG  operation.  The  Nd:YAG 
produced  a  300  mJ  pulse  at  a  10  Hz  repetition  rate  with  a  pulse  FWHM  of  8  ns.  The 
output  pulse  was  periodically  checked  using  various  methods  for  correct  seeder  and 
resonator  alignment  to  ensure  the  output  was  well  seeded  with  a  uniform  transverse 
intensity  profile.  The  first  and  simplest  method  used  bum  paper  to  view  the  beam’s 
intensity  profile.  A  more  precise  method  employed  a  triggered  infrared  camera,  Spiricon 
Pyrocam  I.  The  Nd:YAG  beam  was  highly  attenuated  and  imaged  onto  the  camera.  The 
Nd:YAG  alignment  was  tuned  such  that  the  output  pulse’s  transverse  intensity  profile 
was  Gaussian  .  The  Nd:YAG  laser  did  not  have  to  be  adjusted  often,  and  remained  well 
aligned  for  the  majority  of  the  experiments. 


Figure  3.3.  Seeded  versus  unseeded  NdrYAG  operation.  Notice  the  seeded  pulse 
build-up  time  is  shorter  and  the  pulse  is  smoother. 
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The  seed  laser’s  wavelength  was  set  to  a  specific  rovibrational  transition  energy,  which 
are  tabulated  in  Table  3.1.  The  v(0,3)  pump  transition  energies,  wavelengths,  and  A 
coefficients  for  the  CO  rovibrational  levels  [30]-[31]..  The  Einstein  coefficients  were 
calculated  using  references  [30]  and  [31].  A  Burleigh  WA-1500  Wavemeter  was  used  to 
measure  the  seed  laser’s  wavelength. 


Table  3.1.  The  v(0,3)  pump  transition  energies,  wavelengths,  and  A  coefficients  for 
the  CO  rovibrational  levels  [30]-[31]. 


v  =  0  to  v  =3 

J  lower 

R-Branch  P-Branch 

cmA-1 

nm 

A  (Hz) 

cmA-1 

nm 

A  (Hz) 

6354.1717 

1573.7693 

0.0047 

- 

- 

- 

i 

6357.8064 

1572.8695 

0.0057 

6346.5869 

1575.6500 

0.0137 

2 

6361.3358 

1571.9969 

0.0061 

6342.6372 

1576.6312 

0.0090 

3 

6364.7598 

1571.1512 

0.0065 

6338.5828 

1577.6397 

0.0080 

4 

6368.0782 

1570.3325 

0.0067 

6334.4238 

1578.6756 

0.0075 

5 

6371.2909 

1569.5406 

0.0069 

6330.1604 

1579.7388 

0.0072 

6 

6374.3977 

1568.7757 

0.0071 

6325.7927 

1580.8295 

0.0070 

7 

6377.3984 

1568.0375 

0.0072 

6321.3209 

1581.9478 

0.0068 

8 

6380.2930 

1567.3261 

0.0074 

6316.7452 

1583.0938 

0.0067 

9 

6383.0812 

1566.6415 

0.0075 

6312.0656 

1584.2674 

0.0065 

10 

MSil 

1565.9836 

0.0076 

6307.2824 

1585.4689 

0.0064 

11 

6388.3381 

1565.3523 

0.0078 

6302.3956 

1586.6982 

0.0063 

12 

6390.8065 

1564.7477 

0.0079 

6297.4055 

1587.9555 

0.0062 

13 

6393.1679 

1564.1698 

0.0080 

6292.3121 

1589.2409 

0.0061 

14 

6395.4223 

1563.6184 

0.0082 

6287.1157 

1590.5545 

0.0060 

15 

6397.5694 

1563.0936 

0.0083 

6281.8163 

1591.8963 

0.0059 

16 

6399.6092 

1562.5954 

0.0084 

6276.4142 

1593.2664 

0.0058 

17 

6401.5415 

1562.1237 

0.0086 

6270.9094 

1594.6651 

0.0057 

18 

6403.3662 

1561.6786 

0.0087 

6265.3021 

1596.0922 

0.0056 

19 

6405.0830 

1561.2600 

0.0088 

6259.5924 

1597.5481 

0.0055 

20 

6406.6919 

1560.8679 

0.0089 

6253.7806 

1599.0328 

0.0054 

21 

6408.1927 

1560.5024 

0.0091 

6247.8667 

1600.5463 

0.0054 

22 

BS 

1560.1634 

0.0092 

6241.8509 

1602.0889 

0.0053 

23 

6410.8693 

1559.8509 

0.0093 

6235.7334 

1603.6606 

0.0052 

24 

6412.0449 

1559.5649 

0.0095 

6229.5143 

1605.2616 

0.0051 

25 

6413.1118 

1559.3054 

0.0096 

6223.1937 

1606.8920 

0.0051 

26 

6414.0699 

1559.0725 

0.0097 

6216.7718 

1608.5519 

0.0050 

27 

6414.9190 

1558.8661 

0.0099 

6210.2487 

1610.2415 

0.0049 

28 

6415.6589 

1558.6863 

0.0100 

6203.6246 

1611.9608 

0.0049 

29 

6416.2896 

1558.5331 

0.0102 

6196.8997 

1613.7102 

0.0048 
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The  majority  of  the  Nd:YAG  pulse  energy  was  dumped  to  avoid  damage  to  the 
OPO  optics  using  a  half  wave  plate  and  a  thin  film  polarizer.  The  Nd:YAG  pulse  was 
initially  horizontally  polarized  as  it  exited  the  laser.  Using  a  series  of  polarization 
selective  optics,  Figure  3.2,  the  Nd:YAG  beam  was  selectively  attenuated.  A  maximum 
of  150  mJ  per  pulse  was  sent  into  the  OPO  to  avoid  damage. 

The  attenuated  beam  diameter  was  reduced  through  a  reducing  telescope  to  3  mm, 
which  increased  the  intensity  within  the  nonlinear  crystal.  The  pulse’s  polarization  was 
rotated  to  vertical  orientation  using  another  half  wave  plate  just  before  entering  the  OPO 
resonator. 

CW  Diode  Seed  Laser 

A  tunable  single  mode  laser  diode,  New  Focus,  was  use  to  seed  the  OPO.  The 
diode  was  tunable  between  1.55-1.58  microns  with  an  average  output  power  of  7  mW. 
The  seed  laser  passed  through  two  Faraday  optical  isolators  to  protect  the  diode  from 
unwanted  optical  feedback.  The  single  mode  operation  of  the  New  Focus  diode  was 
verified  using  a  6  GHz  free  spectral  range  etalon  and  an  Indigo  System’s  NIR  camera. 
The  beam  was  imaged  onto  the  NIR  camera  after  passing  through  the  fixed  etalon,  and 
well  formed  rings  were  observed,  indicating  single  mode  operation.  The  New  Focus 
control  box  also  had  a  display  indicating  single  mode  operation. 

OPO  Wavelength  Tuning 

The  OPO  output  pulse’s  wavelength  was  set  by  the  seed  laser’s  wavelength.  The 
seed  laser  forced  the  OPO  to  lase  on  a  single  frequency  due  to  pre-saturation  of  the  gain. 
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OPO  Operation  and  Alignment  Procedure 

The  details  of  the  OPO  will  be  covered  first,  and  then  the  alignment  procedure 
will  be  presented.  The  OPO  cavity  was  a  ring  resonator  design.  See  Figure  3.4  for  a 
detail  layout  of  the  OPO  resonator.  The  output  coupler  was  a  1064  nm  high  transmitter 
and  a  70  %  reflector  at  1550-1580  nm  at  a  22  degree  angle  of  incidence.  The  seed  laser 
was  injected  into  the  cavity  through  the  output  coupler,  and  only  30  %  of  the  seed  laser 
beam  entered  into  the  ring  cavity.  The  back  mirror,  1550-1580  nm  high  reflector,  was 
mounted  on  a  Melles  Griot  piezo-electric  mount.  The  piezo  mount  allowed  for  quick  and 
fine  adjustments  in  the  cavity  length,  and  was  controlled  by  an  Aculight  electronic 
locking  circuit.  The  back  mirror,  high  reflector,  reflects  1550-1580  nm  and  transmits  the 
1064  nm  pump  pulse. 

The  cavity  length  fluctuated  with  room  temperatures,  air  currents,  and  many  other 
outside  variables.  The  locking  circuit  monitors  the  1 .5  pm  light  that  bleeds  through  from 
the  high  reflector,  and  in  turn  adjusts  the  cavity  length  through  the  piezo-electric  mount. 
The  bleed  through  was  maximized  by  the  locking  circuit,  which  occurs  when  the  total 
round  trip  cavity  length  is  equal  to  an  integer  multiple  of  the  wavelength.  The  better  the 
alignment,  the  more  intense  the  bleed  through  becomes. 

Part  of  the  OPO  output  was  split  off  and  directed  through  a  6  GHz  etalon.  The 
rings  were  imaged  on  to  a  CCD  camera  and  observed  on  a  television  monitor.  During 
single  mode  operation  the  rings  were  bright  and  well  spaced.  The  rings  did  not  fluctuate 
in  intensity  or  position  from  pulse  to  pulse.  When  the  OPO  was  misaligned,  not  locked, 
or  not  seeded  the  nonlinear  KTP  crystal  produced  broad  band  1 .5  pm  radiation.  This  was 
evident  by  dimmer,  quickly  shifting  rings  through  the  etalon. 
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OPO  Alignment 

The  OPO  alignment  was  broken  down  into  two  phases.  First  and  most  difficult 
was  aligning  the  seed  laser  through  the  cavity.  Second,  the  Nd:YAG  pulse  must  be 
aligned  to  overlap  the  seed  laser  within  the  nonlinear  KTP  crystal.  A  series  of  carefully 
positioned  apertures  were  used  to  map  out  the  desired  beam  path,  please  reference  Figure 
3.4  for  the  aperture  locations.  The  seed  laser  beam  was  roughly  aligned  through  the 
apertures  by  adjusting  the  two  seed  laser  alignment  optics  (see  Figure  4). 


<y»-.  ft 
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Figure  3.4.  Blow  up  of  OPO  resonator.  Notice  the  aperture  and  camera  location. 

The  apertures  forced  the  beam  to  travel  in  a  straight  and  level  path,  and  the  correct  angle 
through  the  cavity.  A  very  sensitive  InGaAs  near  infrared  camera  (Indigo  Systems  Alpha 
NIR)  was  placed  behind  the  OPO  cavity  high  reflector.  The  seed  laser  injection  optics 
were  adjusted  until  the  seed  laser  enters  the  ring  cavity  and  passes  through  the  aperture  in 
the  center  and  a  circular  diffraction  pattern  is  visible  on  the  InGaAs  camera.  The  seed 
laser  was  now  correctly  aligned  through  the  resonator.  The  output  coupler  was  adjusted 
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until  the  reflected  portion  of  the  seed  beam  travels  through  the  two  widely  spaced 
apertures  (AP#1  and  AP#2  in  Figure  3.4).  The  high  reflector  was  adjusted  such  that  the 
reflected  beam  traveled  through  the  same  widely  spaced  apertures.  This  guaranteed  the 
ring  resonator  seed  beam  is  traveling  through  the  cavity  on  a  straight  and  level  path.  The 
reflections  from  multiple  passes  were  seen  on  the  NIR  camera,  and  the  high  reflector  was 
adjusted  such  that  all  the  reflections  overlapped  one  another.  The  piezo  voltage  was 
manually  adjusted  causing  fluctuations  in  the  cavity  round  trip  length.  The  voltage  was 
tuned  in  such  a  manner  that  the  cavity  goes  through  various  constructive  and  destructive 
interference  lengths.  The  alignment  was  fined  tuned  to  maximize  the  constructive 
interference  visible  on  the  NIR  camera,  see  Figure  3.5. 


Figure  3.5.  Aligned  OPO  bleed  through.  The  typical  beam  profile  display  of  the  bleed 
through  from  the  ring  resonator  for  a  minimized  and  maximized  resonator  length. 


The  Nd:YAG  proved  to  be  a  much  simpler  process.  The  Nd:YAG  pulse  was 
simply  aligned  through  the  ring  cavity  such  that  it  overlapped  the  seed  laser  beam.  The 
Nd:YAG  alignment  was  fine  tuned  once  the  OPO  was  thermally  stabilized  by  tweaking 
the  injection  angle  and  position  to  maximized  output  power.  The  OPO  was  now 
nominally  aligned. 
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The  OPO  Pulse 


The  OPO  typically  produced  horizontally  polarized  8  ns  pulses  with  an  average 
energy  of  15  mJ  at  a  10  Hz  repetition  rate.  The  pulse  energy  was  measured  with  a  Gentec 
ED-200  joule  meter  and  a  Gentec  Dou  reader.  The  data  is  collected  from  the  Dou  meter 
using  the  manufacture’s  Lab  VIEW  software.  The  vast  majority  of  the  Nd:YAG  pulse 
went  unused  and  simply  passed  through  the  cavity  and  was  later  discarded  using  a 
dichroic  optic. 

The  OPO  spatial  beam  profile  was  imaged  using  a  Spiricon  Pyrocam  I  infrared 
camera.  The  spatial  profile  varies  greatly  depending  on  the  OPO  cavity  alignment.  The 
beam  quality  did  not  affect  the  CO  lasing,  however  the  beam  diameter  did  effect  the 
lasing  output  of  the  CO.  This  was  a  direct  effect  of  the  pump  beam  volume  that  was 
incident  on  the  CO  cell.  Since  the  pump  pulse  was  far  from  saturating  the  transition,  the 
same  amount  of  CO  molecules  were  excited  with  a  large  and  small  pump  beam  diameter. 
However,  the  number  density  of  excited  CO  molecules,  directly  related  to  gain,  is 
inversely  related  to  the  pump  volume.  So,  if  the  pump  pulse  was  focused  through  the 
CO  cell,  very  small  beam  radius,  a  large  gain  value  would  be  achieved,  which  will  be 
discussed  further  in  the  results  section.  Figure  3.6  shows  a  typical  spatial  profile  of  the 
OPO  pulse.  CO  lasing  was  observed  under  a  variety  of  pump  pulse  configurations; 
including  collimated  with  a  large  radius  (5  mm),  tightly  collimated  with  a  radius  of 
approximately  2  mm,  and  tightly  focused  in  the  center  of  the  gain  cell.  The  results  and 
effects  will  be  discussed  in  Chapter  4.  Please  see  Figure  3.7  for  a  typical  resonator  and 
pump  pulse  setup. 
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Figure  3.7.  Short  resonator  with  a  double  passed  focused  pump  pulse.  Notice  the 
pump  pulse  slightly  overfills  the  fundamental  mode  volume  of  the  resonator  on  the 
second  pass. 
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Figure  3.8.  OPO  temporal  profile  and  bandwidth.  The  FWHM  of  the  Gaussian  fit  is 
6.0  ns  which  corresponds  to  a  transform  limited  bandwidth  of  166MHz. 

The  OPO  pulse  bandwidth,  Figure  3.8,  is  determined  by  collecting  the  temporal  profile  of 
the  pulse  and  computing  the  Fourier  transform.  The  transform  limited  OPO  pulse 
bandwidth  is  around  150  MHz. 

OPO  Operating  Procedures 

The  OPO  was  a  very  sensitive  optical  system  and  required  great  care  in  setting  its 
operating  parameters  and  conditions.  The  Nd:YAG  was  cycled  on  first  to  a  “flash  only” 
mode  allowing  the  flash  lamps  and  rods  to  come  to  a  thermal  equilibrium.  This  warm  up 
period  was  typically  45  minutes.  While  the  Nd:YAG  was  warming  up  the  rest  of  the 
electronics  were  cycled  on,  oscilloscopes,  lock-in  amplifier,  and  the  piezo  controller. 

The  New  Focus  diode  was  turned  on,  and  the  wavelength  was  adjusted  to  the  correct 


51 


transition  wavelength.  The  diode  also  has  a  30  minute  warm  up  time  according  to  the 
manufacture’s  manual.  Thermal  variations  and  differing  air  currents  from  day  to  day 
caused  the  OPO  to  misalign  over  time.  The  high  reflector  in  the  ring  cavity  was  adjusted 
vertically  and  horizontally  to  re-maximize  the  constructive  interference.  This  re¬ 
alignment  was  performed  daily. 

Once  the  warm  up  time  is  complete  and  the  OPO  is  optimized  the  Nd:YAG  is 
turned  to  Q-Switched  mode  thus  producing  high  energy  pulses.  The  OPO  crystal  is 
allowed  to  thermally  stabilize  for  approximately  45  minutes  before  any  adjustments  are 
made.  Once  the  KTP  crystal  was  stable,  the  OPO  was  fined  tuned  to  maximize  the 
output  pulse  energy  and  pulse  to  pulse  stability. 

3.2  CO  Lasing  Cell 

A  2  inch  diameter  stainless  steel  cell  was  used  for  most  lasing  experiments.  The 
Cell  has  two  CaF2  on  each  end  at  Brewster’s  angle  to  maximize  pump  transmission.  The 
cell  length  was  adjustable  between  29  and  45  inches.  There  are  four  ports  on  the  cell 
used  for  CO  filling  and  fluorescence  measurements. 
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Figure  3.  9.  Carbon  Monoxide  laser  cell  schematic. 
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Another  cell  was  used  for  more  flexible  length  adjustments.  It  was  of  the  same  design  as 
Figure  3.  9  expect  the  cell  diameter  was  1  inch,  thus  allowing  a  standard  piece  of 
stainless  steel  tubing  to  be  inserted  to  adjust  the  length. 

The  cell  was  evacuated  every  test  day  to  about  0.01  Torr  using  a  Welch  vacuum 
roughing  pump.  The  cell  was  then  flushed  with  research  grade  (99.99%  purity)  CO 
several  times  to  ensure  a  low  concentration  of  contaminates.  This  procedure  was  done 
daily  to  minimize  contaminate  gas  effects.  Two  MKS  Baratrons  (10  and  100  Torr)  were 
connected  to  the  cell  to  accurately  measure  pressures  from  1  to  100  Torr.  Both  Baratrons 
were  calibrated  on  a  regular  basis  to  a  NIST  traceable  calibration  stand  by  a  certified 
technician.  The  cell  leak  rate  measured  0.25  Torr/hour,  and  the  cell  was  flushed  through 
an  experimental  series  to  avoid  contamination  effects. 

3.3  Photodetectors 

Multiple  photodiodes  and  photodetectors  were  used  for  the  various  experiments. 
There  were  three  different  wavelength  regions  of  interest,  1064nm,  1550  nm,  and  4.6  pm. 
A  New  Focus  InGaAs  2033  large  area  photo  receiver  was  used  for  all  experiments  that 
require  measurement  of  the  CW  seed  laser,  such  as  the  CO  linewidth  and  absorption 
experiments.  The  New  Focus  photo  detector  has  a  5  mm  x  5  mm  area  and  a  2  ps  rise 
time.  This  detector  was  used  because  it  implemented  an  internal  gain  circuit  that 
produced  an  output  signal  on  the  order  of  volts  which  is  easily  read  using  a  National 
Instruments  DAQ  break  out  board  and  A/D  PCI  card. 

A  ThorLabs  DET1 10  InGaAs  photodetector  was  used  to  accurately  measure  the 
temporal  profile  of  a  short  pulse,  such  as  the  Nd:YAG  pulse.  The  ThorLabs  photo  diode 
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has  a  rise  time  of  5  ns,  active  area  of  0.8mm  x  0.8mm,  and  a  spectral  range  from  700- 
1 800nm.  This  detector  was  mainly  used  as  a  trigger  and  to  measure  the  approximate 
temporal  profile  the  Nd:YAG  pulse. 

Two  Boston  Electronics  PEM-L  photovoltaic  detectors  were  used  as  the  primary 
photodetectors  for  the  CO  laser  pulse  at  4.7  pm.  The  PEM-L  was  sensitive  from  2-12  pm 
with  a  rise  time  of  <0.5  ns.  This  detector  proved  to  be  sensitive  enough  to  measure  the 
CO  laser  pulse,  and  fast  enough  to  for  an  accurate  temporal  profile. 

A  Gentec  QE12SP  joule  meter  was  used  to  measure  the  energy  of  the  4.7  pm 
laser  pulse.  The  QE12SP  came  with  a  manufacture’s  reported  sensitivity  to  7.0  pJ  and  a 
3%  uncertainty  in  power  measurements. 

A  New  Focus  1537  12GHz  photo  receiver  was  used  to  trigger  the  DAQ  system, 
and  obtain  accurate  temporal  profiles  of  the  OPO  pulses.  The  photo  receiver  had  a 
response  time  of  30  ps,  and  an  active  area  radius  of  25  pm. 

3.4  Monochromator  Operation 

A  Digikrom  'A  meter  monochromator  was  used  to  spectrally  separate  the  CO 
lasing  output  into  the  individual  lasing  lines.  A  green  HeNe  laser  (540  nm)  was  used  to 
calibrate  and  align  the  monochromator.  The  monochromator  was  controlled  by  a 
Lab  VIEW  computer  interfaced,  Figure  3.10,  and  its  operation  was  synchronized  with  the 
oscilloscope  data  collection. 
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Figure  3.10.  Screen  shot  of  the  data  acquisition  program.  The  Lab  VIEW  program 
was  used  to  control  all  the  monochromator  settings,  and  collect  scans  from  the 
oscilloscope.  The  program  would  then  compile  the  data  from  multiple  wavelength  scans 
and  save  the  data  in  Excel  format. 


3.5  Data  Acquisition:  Oscilloscope 

A  Tektronix  TDS  680B  oscilloscope  was  used  for  the  data  collection  from  the 
various  photo  detectors.  The  TDS  680B  had  a  5  Gs/sec  maximum  sampling  rate  and  a  1 
GHz  bandwidth. 
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Figure  3.11.  Experimental  layout  drawn  to  scale  to  clarify  timing  issues  with  the 
trigger  pulse  and  CO  laser  pulse. 

The  New  Focus  photo  detector  triggered  the  oscilloscope  from  the  scatter  off  a  steering 
mirror,  see  Figure  3.11.  A  1.5  pm  bandpass  filter  ensured  the  photo  detector  triggered 
off  the  1 ,5pm  pulse  of  interest,  and  not  scattered  radiation  from  the  OPO. 

Timing  considerations  were  addressed  early  on  in  the  experiment  by 
standardizing  the  trigger’s  physical  location,  cable  lengths,  and  the  CO  detector 
locations.  A  6  ns  delay  was  experimentally  measured  between  the  trigger  pulse  and  the 
pump  pulse’s  entrance  in  the  resonator.  This  delay  includes  both  cable  and  light 
propagation  delays.  An  1  Ins  delay  was  calculated  between  the  CO  pulse  exiting  the 
resonator  and  reaching  the  PD  #1,  and  a  13  ns  delay  to  PD#2.  The  delay  using  a  long 
resonator  was  6ns  and  7ns  respectively.  All  plots  presented  are  corrected  for  the  timing 
delays  unless  otherwise  noted. 

The  CO  pulses  were  typically  averaged  over  5  waveforms  to  minimize  shot  to 
shot  noise,  and  the  data  were  sent  to  a  computer  using  Lab  VIEW  software  and  a  GPIB 
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interface.  The  data  acquisition  software  displayed  the  plots  on  the  screen  allowing  the 
user  to  view  the  collected  plots  before  saving  them  to  an  Excel  file. 

The  shot  to  shot  variations  depend  directly  upon  the  pulse  stability  of  the  OPO. 
The  OPO  was  not  very  stable  from  shot  to  shot,  but  it  did  produce  a  stable  output  over  5 
to  10  pulses,  Figure  3.12.  A  Gentec  Duo  power  meter  provided  the  average  energy  over 
a  user  defined  number  of  pulses.  So,  the  average  pump  pulse  energy  is  reported  along 
with  the  average  CO  laser  pulse  data. 


Averaging  Effects 
Arbitrary  CO  Laser  Pulse 


Time  (ns) 


Figure  3.12.  Averaging  effects.  The  pink  scan  is  an  average  of  5  scans,  and  the  blue 
trace  contains  no  averaging.  Average  effects  were  considered  through  out  the  experiment 
and  were  considered  negligible. 
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CHAPTER  FOUR:  RESULTS  AND  DISCUSSION 


This  chapter  gives  the  details  and  results  of  the  CO  lasing  experiments,  and  the 
theoretical  modeling. 

4.1  Lineshape  Analysis  and  Pressure  Broadening  Coefficients 

The  2nd  overtone  absorption  lineshapes  of  CO  were  measured  using  the  tunable 
diode  laser  over  range  of  pressures,  from  0-100  Torr.  The  experimental  data  agreed 
within  15%  of  the  published  data  by  Henningsen  [22].  The  difference  can  be  explaining 
by  systematic  errors  in  the  experiment,  such  as  using  a  non  optimized  Fabry  Perot 
interferometer  to  calculate  the  frequency  scan  range  of  the  diode.  The  interferometer  was 
originally  coated  for  wavelengths  ranging  from  1.1  pm  to  1.3  pm,  and  the  finesse 
suffered  greatly  at  1.5  pm.  This  created  a  large  uncertainty  in  the  frequency  calibration 
of  the  laser  scan.  A  simpler  less  stringent  method  was  used  to  fill  the  CO  cell  than  that  of 
Henningsen,  possibly  resulting  in  a  higher  contaminate  ratio.  However,  regardless  of  the 
experimental  error  a  good  agreement  between  published  results  was  obtained, 

Table  2.4.  See  Figure  4.1  for  a  sample  of  the  fitted  lineshape  function  at  10  and  100  Torr 
of  CO. 
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Figure  4.1.  The  lineshapes  of  v(0,3)R(8)  transition  at  10  and  100  Torr  of  CO  with 
pressure  broadened  linewidths  of  260MHz  and  26MHz  respectively. 

4.2  CO  Lasing  Results  Short  Gain 

The  results  of  a  short  gain  medium,  33cm,  will  be  covered  in  this  section.  The 
effect  of  CO  pressure,  output  coupler  transmission,  and  resonator  length  were  measured 
for  the  v(0,3)  R(6)  pump  transition,  as  shown  in  Figure  4.2. 
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Figure  4.2.  Energy  level  diagram  of  pumped  transition. 

A  stable  resonator  with  a  45  cm  gain  length,  a  high  reflector  with  a  5  m  radius  of 
curvature,  and  a  flat  output  coupler  with  reflectivities  from  10-95%  were  used  for  this 
series  of  experiments.  The  resonator  parameters  for  the  short  cavity  can  be  found  in 
Figure  4.3. 
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Figure  4.3.  Fundamental  resonator  mode  and  pump  pulse.  The  blue  lines  represent 
the  resonator  optics,  and  the  solid  red  line  represents  the  double  pass  optic. 


Notice  in  Figure  4.3  the  pump  beam  will  excite  CO  outside  of  the  fundamental 
mode  volume  thus  inducing  higher  order  transverse  mode  lasing.  The  multi-mode  lasing 
caused  an  extreme  sensitivity  on  the  photo  detector  location,  since  the  detector  would 
pick  up  different  mode  beats  and  entirely  different  modes  based  on  its  location.  This 
resulted  in  a  strange  temporal  profile  of  the  CO  laser  pulse.  This  phenomenon  was 
observed  in  the  data  sets,  and  can  be  explained  by  transverse  mode  beats. 

Two  data  sets  were  collected  along  with  the  4.7  pm  beam  profile.  In  both  data 
sets  the  laser  ran  in  multimode  operation,  but  the  resonator  was  misaligned  in  the  second 
data  set  to  allow  for  higher  order  mode  oscillation,  as  shown  in  Figure  4.4. 
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Figure  4.4.  Multi-transverse  mode  lasing.  The  red  line,  MM2,  corresponds  to  the 
bottom  beam  profile,  and  the  blue  line  corresponds  to  the  top  beam  profile. 

The  red  line,  in  Figure  4.4,  has  a  secondary  beat  frequency  that  is  absent  from  the  blue 
line,  due  to  more  transverse  modes  running.  A  qualitative  estimate  of  the  modes  running 
in  each  case  can  also  be  inferred  from  the  beam  profiles  in  Figure  4.4.  The  fast  Fourier 
transform,  FFT,  of  MM1  and  MM2  (Figure  4.4),  was  taken  to  determine  the  beat 
frequencies  more  accurately.  The  beats  line  up  with  the  transverse  mode  frequency 
spacing  between  the  fundamental  mode  and  the  TEMoi  and  TEMn  that  were  calculated 
using  equations  2.62  and  2.63. 
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Figure  4.5.  FFT  of  the  multi-mode  CO  laser  pulse.  The  beat  frequencies  occur  at  the 
two  transverse  mode  beats. 

The  CO  laser  was  forced  to  run  in  single  mode  operation  simply  by  placing  an 
adjustable  aperture  within  the  resonator  and  which  filtered  out  all  higher  order  modes. 
Figure  4.  6  shows  a  plot  of  a  unapertured  resonator  and  an  apertured  resonator.  The 
oscillations  are  no  longer  present  under  the  single  mode  operation  since  only  one  mode 
exists  and  there  are  no  beats.  The  mode  beating  effects  are  evident  in  most  of  the  short 
gain  data  sets. 
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Figure  4.  6.  Multimode  and  single  mode  operation.  The  laser  was  forced  to  run  in  the 
fundamental  mode  by  aperturing  out  all  other  modes  in  the  resonator. 

Data  was  collected  using  two  pump  configurations,  focused  pump  and  unfocused  pump. 
The  focused  pump’s  beam  profile  is  displayed  in  Figure  4.3,  and  the  unfocused  pump 
maintained  an  average  beam  diameter  of  6  mm  throughout  the  gain  region.  The  gain  was 
calculated  based  on  the  percentage  of  the  pump  absorption  and  the  mode  volume  of  the 
pump.  The  gain  volume  was  calculated  as  the  product  of  the  gain  length  and  the  cross- 
sectional  area  of  the  pump  mode.  The  absorption  of  the  pump  pulse  proved  to  be  very 
difficult  to  measure  due  to  fluctuations  in  pump  pulse  intensities  and  the  very  small 
absorption  cross-section  of  the  2nd  overtone.  The  absorption  was  measured 
experimentally  by  simultaneously  using  two  Joule  meters.  The  first  energy  meter 
measured  a  portion  of  the  pump  pulse  not  transversing  the  CO  cell,  and  the  second 
energy  meter  measured  the  power  transmitted  through  the  cell.  The  ratio  between  these 
two  powers  any  noise  from  pump  pulse  fluctuations.  The  CO  cell  was  then  filled  with 
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various  pressures  of  CO  and  the  ratios  were  compared  to  that  of  the  empty  cell.  The 
results  are  shown  in  Figure  4.7. 


Figure  4.7.  Pump  absorption  versus  CO  pressure  through  a  1.82m  static  CO  cell. 

The  absorption  curves  above  are  for  the  v(0,3)R(6)  transition  at  1567.77nm. 

The  gains  at  30  Torr  of  CO  for  the  focused  and  unfocused  pump  profile  was  estimated  to 
be  around  150  %  cm'1  and  30  %  cm'1  respectively,  with  out  taking  into  account  any  loss 
from  R-R.  The  huge  difference  in  the  gains  was  due  to  the  squared  dependence  of  gain 
on  the  pump  radius.  Under  the  focused  pump  conditions  the  high  gain  results  in  an  ASE 
pulse  allowing  the  CO  to  lase  without  any  resonator  optics.  A  typical  set  of  spectrally 
resolved  plots  are  shown  in  Figure  4.8.  Some  salient  observations  were  made,  such  as 
the  mode  build  up  times  between  the  focused  and  unfocused  pump  pulse.  The  focused 
pump  pulse  has  a  mode  build  up  time  of  5ns  and  was  coincident  with  the  pump  pulse. 
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This  was  evident  of  the  fact  that  the  CO  was  super  radiant  under  the  same  conditions. 
The  unfocused  pump  case  had  a  mode  build  up  time  on  the  order  of  15ns,  and  was  not 
super  radiant,  Figure  4.8.  The  value  agrees  with  the  results  from  the  analytical 
representation  of  the  mode  build  up  time  [5] 

r  =  — In  {ALNVS)  (4.1) 

cLa 

where  d  is  the  resonator  length,  L  is  the  gain  length,  a  is  the  small  signal  gain,  A  is  the 
mode  area,  and  8  is  used  to  represent  the  fractional  change  in  gain  (£=0.1).  Equation 
4.1  was  used  to  determine  the  mode  build  up  time  based  off  of  the  calculated  gain,  and 
also  to  calculate  the  gain  from  the  measured  mode  build  up  time.  The  numerically 
estimated  mode  build  up  times  for  30  Torr  of  CO  the  unfocused  pump  is  16.2  ns.  This 
number  agrees  very  well  with  the  experimentally  determined  value  of  14  ns  (unfocused). 
Equation  4. 1  was  solved  for  gain  and  the  upper  state  number  density  simultaneously 
using  the  measured  mode  build  up  times.  Gain  calculated  from  equation  4.1  for  the 
unfocused  pump  was  10  %  cm'1.  This  compares  well  with  the  gain  calculated  based  on 
the  pump  absorption  and  pump  beam  area,  7  %  cm'1. 
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Figure  4.8.  A  typical  plot  of  the  spectrally  resolved  CO  laser  pulse.  The  oscillations 
are  described  as  transverse  mode  beats.  The  lasing  only  occurred  on  the  2nd  overtone 
band.  This  data  represents  an  unfocused  pump  pulse  with  a  gain  on  the  order  of  30  % 
cm'1. 

Cascade  lasing,  lasing  from  the  v(3,2)  and  v(2,l)  levels,  was  not  seen  under  any  resonator 
configuration  or  CO  pressure  for  the  unfocused  pump  pulse.  This  behavior  was  simply  a 
result  of  the  much  lower  gain.  Population  inversion  initially  starts  on  the  2nd  overtone 
and  it  lases  to  the  1st  overtone  level.  Inversion  on  the  v(2,l)  line  was  created  by  the 
radiative  stimulated  emission  from  v(3,2).  However,  the  initial  inversion  is  split  by  the  P 
and  R  branch  lasing  lines  to  two  separate  rotational  levels  in  the  1st  overtone  band.  This 
result  effectively  implies  the  inversion  on  the  v(2,l)  transition  is  at  best  %  of  the  initial 
inversion.  The  smaller  inversion  of  the  1st  overtone  level  causes  the  gain  to  be  below  the 
threshold,  and  stimulated  emission  never  begins. 
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Figure  4.9.  Typical  plot  of  the  spectrally  resolved  CO  pulse  under  the  focused  pump 
transition.  The  gain  is  much  higher  under  these  conditions,  and  cascade  lasing  is  evident 
(short  cell). 


The  short  gain  cell  condition  was  not  investigated  in  great  detail  because  it  was  far  from 
an  optimized  condition.  A  very  small  fraction  of  the  pump  pulse  was  absorbed  in  the 
small  gain  length.  It  was  used  to  experimentally  show  that  under  a  long  resonator  and 
long  gain  cell  conditions,  on  the  order  of  2  meters,  the  CO  laser  self  mode-locks.  Since 
the  FSR  of  the  short  resonator  is  on  the  order  of  270  MHz  and  the  lasing  transition’s 
FWHM  is  also  on  the  order  of  270  MHz  only  one  axial  mode  would  be  present.  This 
would  eliminate  the  CO  laser’s  mode-locking  capability  due  to  the  lack  of  multiple 
modes. 
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4.3  CO  Lasing  Results  Long  Gain 

A  long  gain  cell  of  2.0  m  was  used  for  the  vast  majority  of  the  data  collection. 
Under  the  long  gain  cell  a  variety  of  conditions  were  studied.  An  unexpected  result 
occurred  with  the  long  gain  cell  and  a  focused  pump  beam  geometry,  the  CO  laser  self 
mode-locked. 

The  CO  was  super  radiant  under  these  conditions  and  did  not  require  a  resonator 
to  emit  a  semi  coherent  pulse.  The  most  interesting  part  was  that  the  super  radiant  pulse 
acted  as  an  injection  seeder  and  forced  the  CO  laser  to  mode-lock,  Figure  4.10. 


Figure  4.10.  Mode-Locked  pulse  train. 

The  v(3,2)P(8)  lasing  line  has  a  FHWM  of  400  MHz  at  a  CO  pressure  of  75  Torr, 
and  the  resonator  had  a  full  spectral  range  of  74  MHz.  This  implies  approximately  5 
axial  modes  are  resonating  within  the  cavity  if  a  highly  reflective  output  coupler,  95% 
reflective,  is  used.  The  pulse  widths  are  related  to  the  number  of  modes  locked  together 
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The  mode-locking  operation  was  evident  only  under  the  focused  pump  conditions, 
higher  CO  pressures,  and  a  high  reflectivity  output  coupler.  The  high  output  coupler’s 
reflectivity  lowers  the  threshold  gain  level,  and  allows  for  more  axial  modes  to  break 
threshold.  The  focused  pump  condition  increased  the  gain  which  also  allows  for  more 
axial  modes  to  break  threshold.  The  higher  CO  pressures  broadened  the  CO  lasing  line 
spectral  width,  again  allowing  more  axial  modes  to  fit  within  the  gain  profile. 

The  next  series  of  experiments  were  designed  to  maximize  the  CO  laser  power 
output,  and  efficiency.  A  long  gain  cell,  1 82cm,  was  used  along  with  a  long  resonator, 
233  cm.  Specifications  for  this  series  of  experiments  can  be  found  in  Table  4.  1 . 


Table  4. 1.  Experimental  specifications  for  lasing  experiments. 


CO  Laser  Specifications 

182  cm 

Resonator  Length 

233  cm 

Cavity  Round  Trip 

15.5ns 

Cavity  FSR 

64.4MHz 

MaxR 

98%  R 

Radius  of  Curvature 

500  cm 

Pump  Transition 

v(0,3)R(6) 

Pump  Wavelength 

1568.77nm 

The  pump  pulse  was  double  passed  through  the  static  CO  cell.  Various  CO  pressures, 
pump  beam  profiles,  and  output  coupler  reflectivities  were  explored.  The  pump  beam 
geometry  used  to  produce  the  maximum  power  was  a  tightly  collimated  pump  pulse  with 
an  average  beam  waist  of  1.89  mm  through  the  gain  region  including  the  double  pass. 
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Figure  4.12.  Long  cavity  resonator  fundamental  mode  volume.  It  is  easily  seen  that 
the  pump  beam  creates  gain  only  within  the  cavity’s  fundamental  mode  volume.  This 
implies  the  CO  laser  is  running  in  single  mode  operation. 

ASE  was  present  under  the  optimized  pump  profile;  however  it  was  not  as 
dominating  as  under  the  focused  conditions.  Mode-locking  was  also  less  prevalent  in  the 
optimized  geometry.  It  appeared  as  though  the  CO  laser  was  running  only  in  the 
fundamental  transverse  mode,  since  no  transverse  mode  beat  frequency  appeared  in  the 
data  and  the  pump  was  mode  matched  to  the  fundamental  mode,  Figure  4.12.  However, 
this  was  not  experimentally  confirmed. 


72 


Figure  4.13.  CO  laser  pulse  out  with  a  70%  R  output  coupler  under  various  CO 
pressures. 

The  CO  pressure  had  a  large  effect  on  the  temporal  behavior  of  the  CO  laser  pulse, 

Figure  4.13.  At  lower  pressures  the  CO  lased  for  a  much  longer  period  of  time,  but  with 
lower  peak  intensity.  This  observed  behavior  can  be  explained  by  the  RR,  W,  and  VT 
effects  that  start  to  come  into  consideration  at  higher  CO  pressures.  As  the  pressures 
become  higher  the  kinetic  effects’  time  scales  decreases  to  that  of  the  lasing  time  scales. 
All  the  above  kinetic  effects  decrease  the  population  inversion  with  time,  thus  lowering 
the  gain.  The  initial  spike  that  overlaps  with  the  pump  pulse  is  the  ASE  pulse,  which  is 
larger  at  lower  pressures.  The  gain  calculated  from  the  absorption  and  mode  build  up 
time  was  estimated  to  be  22  ±  1  %  cm'1.  The  optimized  pressure  was  found 
experimentally  to  be  between  40  and  50  Torr  of  CO  for  all  the  output  coupler 
reflectivities.  The  CO  laser  pulse  was  spectrally  resolved  through  a  monochromator,  and 
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various  lasing  lines  were  observed.  The  CO  did  not  lase  on  the  v(2,l)  transition  at  higher 
pressures,  above  20  Torr.  This  was  due  to  the  RR  and  W  decreasing  the  gain  on  the 
specific  transition.  Spectrally  resolved  lasing  transitions  with  a  70%  reflective  output 
coupler  can  be  seen  in  Figures  4.14-15  at  various  CO  pressures.  At  pressures  above  20 
Torr  no  cascade  lasing  was  present.  Please  see  Appendix  A  for  the  plots  of  all  the  data 
collected. 


Time  (ns) 


Figure  4. 14.  Figure  4.14.  Spectrally  resolved  lasing  output  with  a  70%  reflective 
output  coupler,  and  10  Torr  of  CO.  There  are  two  lasing  lines  from  the  v(3,2)  both  the 
P  and  R  lines,  and  two  v(2,l)  only  the  P  lines. 


74 


Intensity  (A.U.) 


Figure  4.15.  Spectrally  resolved  lasing  output  with  a  70%  reflective  output  coupler, 
and  20  Torr  of  CO. 

An  intensity  saturation  curve  was  collected  by  measuring  the  CO  laser  pulse 
energy  versus  output  coupler  transmission,  the  results  can  be  seen  in  Figure  4.16.  The 
Rigrod  analysis  was  applied  to  the  collected  Isat  curves,  however  the  results  have  very 
limited  accuracy  due  to  several  reasons.  First,  the  Rigrod  analysis  assumes  continuous 
wave  operation  which  is  far  from  the  experimental  setup.  Since,  the  pulse  width  is  on  the 
order  of  6ns  and  the  cavity  round  trip  time  is  around  15ns  the  pump  pulse  does  not  engulf 
the  entire  cavity  at  one  time,  rather  it  rattles  back  and  forth.  This  causes  the  gain  change 
spatially  and  temporally.  Also,  the  Rigrod  analysis  does  not  take  into  account  ASE, 
which  also  affects  the  accuracy  of  the  prediction. 
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Output  Coupler  Transmission  (%) 


Figure  4.16.  Experimental  output  energy  versus  output  coupler  transmission. 

The  effect  of  pressure  on  lasing  efficiency  was  also  studied  experimentally  and 
theoretically.  A  comparison  between  the  numerical  model  and  experimental  results  is 
presented  in  Figure  4.18.  The  numerical  model  agrees  well  with  the  experimental  data. 
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Figure  4.17.  Power  output  versus  pressure  with  a  40%  reflective  output  coupler. 

4.4  Discussion 

Overall  the  experiments  were  a  success  and  optically  pumped  cascade  lasing 
around  4.7  um  was  demonstrated  on  the  2nd  overtone  of  carbon  monoxide.  The  data 
agreed  well  with  the  numerical  model  results,  Figure  4.17,  and  the  max  optical  efficiency 
measured  was  12  ±  1%.  The  low  optical  efficiency  can  be  partially  explained  by  to 
cavity  losses  and  older  optics  used  in  the  resonator.  Energy  efficiency,  defined  as  the 
ratio  of  pump  energy  input  to  the  laser  energy  extracted  from  the  resonator,  was 
experimentally  around  2  ±1%.  The  low  energy  efficiency  can  be  attributed  to  low  pump 
absorption  (a  maximum  of  25%  of  the  pump  was  absorbed).  The  numerical  model 
suggests  overall  energy  efficiency  at  room  temperature  using  a  double  pass  2  m  cell  can 
approach  8  %  with  a  high  Q-cavity,  with  an  optical  efficiency  of  33%.  As  the  pump 
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absorption  approaches  100  %,  both  the  energy  efficiency  and  optical  efficiencies  improve 
drastically  to  32%. 

The  model  lacked  the  ability  to  accurately  describe  the  temporal  output  of  the  CO 
laser  since  it  did  not  incorporate  the  resonator  modes  into  the  calculations.  However,  it 
proved  useful  in  determining  overall  lasing  trends.  The  model  also  lacked  an  ASE  term, 
which  also  accounted  for  some  the  discrepancies  between  predictions  and  measurements. 

Future  experimental  investigation  should  include  cooling  the  CO  cell  to  100K  in 
order  to  increase  the  absorption  efficiency.  The  CO  should  also  be  heated,  which  has 
known  draw  backs.  However,  the  higher  v(0,3)  rotational  levels  have  energy  gaps  closer 
to  the  1550  nm  wavelength  of  Erbium  fiber  lasers.  Another  improvement  would  be  to 
fine  tune  the  OPO  and  amplify  the  1567  nm  pump  pulse.  These  steps  would  increase  the 
photon  flux  through  the  CO  cell,  and  improve  upon  the  power  output.  The  numerical 
model  could  be  improved  with  the  addition  an  ASE  term. 

The  data  did  not  show  much  benefit  in  cascade  lasing,  since  the  highest  powers 
measured  were  without  cascade  lasing.  This  leads  to  the  question  of  why  bother  with  the 
small  absorption  cross-sections  of  2nd  overtone  pumping  and  the  absorption  problems  it 
poses.  Instead,  a  2.3  pm  GaSb  diode  laser  could  be  used  as  a  1st  overtone  optically 
pumped  source.  The  absorption  cross  section  at  2.3  pm  is  over  100  times  larger.  The  2.3 
pm  diode  laser  is  not  at  the  same  developmental  stage  as  the  Erbium  fiber  lasers,  and 
lacks  the  high  power  to  produce  a  significant  4.7  pm  laser. 
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CHAPTER  FIVE:  FUTURE  STUDY 


The  method  used  to  create  an  optically  pumped  CO  laser  is  inherently  unscalable. 
Since,  OPO’s  do  not  scale  linearly  and  are  non-trivial  devices  to  package  in  a  turn  key 
system  the  pump  source  is  currently  the  limiting  factor  on  optically  pump  CO  overtone 
laser’s  efficiency  and  power  level.  In  order  to  increase  the  CO  laser’s  output  power  a 
more  powerful  pump  source  must  be  explored.  Currently,  high  power  fiber  lasers  are 
developing  dramitically,  and  have  already  surpassed  the  lkW  average  power  per  fiber 
level.  The  idea  of  using  fiber  lasers  to  optically  pump  CO  will  be  briefly  explored. 

Fiber  lasers  currently  pose  two  difficulties  in  optically  pumping  CO.  First,  they 
operate  at  extremely  high  repetition  rates  which  can  be  considered  continuous  wave,  CW, 
operation.  Second,  they  typically  have  a  very  large  bandwidth  range.  The  large 
bandwidth  decreases  the  absorption  efficiency.  Research  is  currently  underway  to 
narrow  band  a  high  power  fiber  laser  around  the  v(0,3)  transition  wavelength,  1567nm. 

CW  operation  is  a  problem  since  the  intensity  drops  drastically  compared  to  pulse 
operation.  This  behavior  hints  the  power  levels  required  to  adequately  pump  the  2nd 
overtone  of  CO  with  fiber  lasers  may  be  unreasonably  high.  The  numerical  model 
described  in  Appendix  C  was  used  to  run  a  numerical  study  on  CW  pumping  of  CO. 

A  long  trapezoidal  pulse  was  used  to  approximate  a  CW  pump  source.  The  pulse 
length  was  set  to  2000  ns,  which  was  long  enough  for  the  CO  laser  to  reach  steady  state, 
as  shown  Figure  5.1. 
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Figure  5.1.  CO  laser  output  with  a  25kW  CW  pump  source.  The  green  trace  is  the 
scaled  pump,  and  the  orange  trace  is  the  total  laser  output.  The  blur  of  lines  are  high 
rotational  level  lasing  lines  due  to  a  partial  population  inversion. 


The  model  outputs  the  absorption  efficiency  and  the  energy  efficiency.  Absorption 
efficiency  is  defined  as  the  percentage  of  the  pump  beam  absorbed,  and  the  energy 
efficiency  is  the  fraction  of  pump  energy  over  CO  energy  extracted  from  the  resonator. 
Table  5.1  has  a  list  of  the  model  parameters  used  to  pump  the  v(0,3)  R(2)  line. 


Table  5.1.  Model  input  conditions  for  CW  pumping  runs. 


MODEL  PARAMETERS 

CO  Pressure 

30  Torr 

100  cm 

110  cm 

90%  R 

0.30  cm 

100  MHz 

50 
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The  first  parameter  explored  was  the  effect  on  CW  pump  power  versus  CO  laser 
efficiency.  Model  runs  were  done  at  two  temperatures,  298K  and  100K.  At  lower 
temperatures  a  smaller  Boltzmann  distribution  is  obtained,  see  Figure  5.2.  These 
conditions  raise  the  efficiency  by  increasing  the  amount  of  the  pump  that  can  be  absorbed 
before  the  pump  transition  is  saturated,  and  the  lower  rotational  states  have  a  much  larger 
cross  section. 


Figure  5.2.  Boltzmann  rotational  distribution  for  the  ground  state  of  CO.  Notice  the 
number  of  molecules  in  J  =  4  at  100  K  is  approximately  double  that  of  J  =  6  at  298  K. 

The  increased  absorption  cross  section  allows  a  higher  percentage  of  the  pump  beam  to 
be  absorbed.  In  other  words  it  reduces  the  effective  gain  length  needed  to  obtain  100% 
absorption. 
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The  pump  beam  is  97%  at  298  K  and  99%  at  100  K  absorbed  in  the  analyses  represented 
in  Figure  5.3.  The  pump  beam  saturates  the  transition  around  7.5  kW  at  298  K.  The 
pump  does  not  saturate  the  transition  at  100  K  even  at  powers  over  50  kW,  as  shown  in 
Figure  5.4. 
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Figure  5.4.  Absorption  versus  pump  power.  This  shows  the  effect  of  saturation  on  the 
pump  beam  absorption.  Running  at  cooler  temperatures  allows  for  more  of  the  pump  to 
be  absorbed,  thus  increasing  the  efficiencies. 


Figure  5.  5.  Absorption  versus  path  length  with  30  Torr  of  CO.  At  100K  99%  of  the 

beam  is  absorbed  within  10  meters,  at  298K  50  meters  is  required. 
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An  intensity  saturation  curve,  Figure  5.  6,  was  modeled  under  the  conditions  in 
Table  5.1.  Results  indicate  that  the  gain  is  still  high  in  the  partial  inversion  since  the 
model  predicts  lasing  at  very  high  output  coupler  transmissions. 


Figure  5.  6.  The  Isat  curve  for  30  Torr  of  CO  at  298K  and  100K. 

The  numerical  results  show  promise  in  the  ability  to  CW  pump  CO,  with  very 
high  efficiencies.  The  model  indicates  the  small  absorption  on  the  2nd  overtone  can  be 
overcome  by  cooling  the  CO  to  100K  reducing  the  length  required  from  50  m  to  10  m  for 
99%  pump  absorption.  The  overall  energy  efficiency  approaches  50%  for  CW  powers 
over  12kW. 


APPENDIX  A:  DATA  PLOTS 
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This  appendix  contains  all  of  the  data  collected  for  the  long  gain  and  resonator 
setup.  The  specifications  listed  in  Table  A.l  are  valid  for  all  plots  in  appendix  A. 

Table  A.  1.  Specifications  for  all  data  presented  in  Appendix  A. 


CO  Laser  Specifications 

Gain  Length 

182  cm 

Resonator  Length 

233  cm 

15.5ns 

Cavity  FSR 

64.4MHz 

MaxR 

98%  R 

Radius  of  Curvature 

500  cm 

Pump  Transition 

Pump  Wavelength 

1568.77nm 
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Figure  A.2.  70%  reflective  output  coupler  and  100  -  5  Torr  of  CO, 
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Figure  A.3.  60%  reflective  output  coupler  and  100-5  Torr  of  CO. 
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Figure  A.4. 40%  reflective  output  coupler  and  100  -  5  Torr  of  CO. 
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Intensity  (A.U.) 
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Figure  A.  6.  ASE  emission  plots  for  various  pressures. 
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APPENDIX  B:  LINE  WIDTH  EXPERIMENTS 


Linewidth  experiments  used  the  OPO’s  diode  seed  laser  as  a  tunable  probe  over 
the  v(0,3)  transition.  The  scanning  laser  output  was  controlled  and  data  collected  using  a 
computer  employing  Lab  VIEW  software  and  National  Instruments  DAQ  hardware.  The 
computer  controlled  the  output  wavelengths,  scan  range,  and  the  scan  rate.  The  central 
output  wavelength  was  set  by  coarsely  tuning  an  intracavity  diffraction  grating  of  the 
diode,  and  the  fast  tight  scans  were  accomplished  by  a  piezo  electric  motor  that  fine 
tuned  the  wavelength  out.  The  piezo  was  driven  by  a  voltage  ramp  generated  by  the 
Lab  VIEW  control  software.  The  scan  rate  was  adjusted  by  recycling  the  voltage  rise 
across  the  diode  at  rates  of  up  to  50  Hz.  Resolution  on  the  frequency  control  was  500 
samples  per  scan. 

The  diode  laser  beam  was  split  two  ways.  The  input  measurements  where  data 
collected  on  a  reference  cell  was  used  to  determine  if  the  correct  wavelength  range  was 
being  scanned,  and  Fabry-Perot  Interferometer  to  correctly  calibrate  the  x-axis  to  units  of 
frequency. 

A  Fabry-Perot  Interferometer  was  used  to  correctly  calibrate  the  x-axis  to  one  of 
frequency  versus  samples.  Data  was  read  into  the  system  based  only  on  a  per  sample 
calibration.  For  the  results  of  the  absorption  measurements  to  have  any  meaning,  the  x- 
axis  must  be  calibrated  to  that  of  frequency.  In  order  to  do  so,  it  was  necessary  to  read 
the  output  of  the  laser  diode  on  a  Fabry-Perot  Interferometer.  Since  the  scan  range  was 
approximately  2.4  GHz,  a  300  MHz  Fabry-Perot  Interferometer  was  used  to  ensure 
enough  peaks  could  be  measured  for  adequate  calibration.  Using  the  known  value  of 
300MHz  for  a  full  spectral  range  peak  and  comparing  that  value  to  the  number  of 
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samples  between  peaks,  a  valid  value  of  frequency  could  be  measured  on  the  data 
acquisition  readout. 

Lab  VIEW  software  then  collected  the  absorption  with  an  accurate  frequency  axis, 
and  saved  the  data  to  a  text  file.  This  data  was  then  fit  to  a  Voigt  function  using  Orgin’s 
software  package.  In  order  to  get  accurate  Gaussian  and  Lorentzian  widths  the 
relationships  temperature  and  pressure  must  be  maintained.  The  Gaussian  linewidth  is 
related  only  to  the  temperature  of  the  gas  by 


where  wG  is  the  Gaussian  (Doppler)  width,  vQ  is  the  line  center  frequency,  T  is  the 

temperature,  and  M  is  the  reduced  mass.  The  Lorentzian  width  is  related  to  both  the 
temperature  and  pressure  of  the  gas  by 

(2) 

where  wL  is  the  Lorentzian  width,  y  is  the  pressure  broadening  coefficient  [MHz/Torr] 
of  species  i,  and  P  is  the  pressure  of  the  species  i.  When  fitting  with  a  Voigt  function  the 
resulting  Gaussian  and  Lorentzian  width  must  be  consistent  [33].  This  simply  means  the 
temperature  corresponding  to  both  widths  must  agree. 

Using  the  above  procedure  to  fit  absorption  curves  at  various  pressures  of  CO  the 
pressure  broadening  coefficient  can  be  experimentally  measured. 
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APPENDIX  C:  Numerical  Model  Description 

The  rate  equation  model  was  developed  using  existing  code  to  model  the  optically 
pumped  CO  laser.  The  model  included  terms  for  rotation  relaxation,  vibrational- 
vibrational,  vibrational-rotational,  and  radiative  processes.  The  model  did  not  originally 
contain  any  spatial  dependence,  longitudinal  or  radial,  and  simply  assumed  the  entire 
gain  region  changed  only  in  time.  The  model  employed  a  robust  ordinary  differential 
equation  solver  to  solve  the  120  coupled  differential  equations.  This  section  heavily 
references  the  model  developed  by  Rafferty  et.  al  [32]. 

The  rate  equations  used  can  be  written  as  one  equation  with  rate  constants  for 
each  individual  rovibrational  level.  This  results  in  120  (30  rotational  levels  and  4 
vibrational  levels)  coupled  differential  equations.  Rate  equations  are  based  on  of  Chapter 
2  theoretical  discussion  of  the  radiative  processes  and  their  effects  on  the  population 
inversions. 

The  model  will  be  developed  here  under  two  processes,  radiative  and  non- 
radiative.  Radiative  processes  include  spontaneous  and  stimulated  emission,  and  the 
non-radiative  processes  include  the  above  mentioned  kinetics.  The  time  rate  of  change  of 
the  number  of  photons  in  the  cavity  is  dependent  upon  the  gain  and  is  given  by 

~r=~r^a~a^Np  0) 

at  a 

where  Np  is  the  number  of  photons,  a  is  the  gain  between  the  upper  and  lower  levels, 
ath  is  the  threshold  gain,  d  is  the  resonator  length,  and  Lg  is  the  gain  length.  It  is 
important  to  notice  that  if  Np  =  0  at  t  =  0,  then  for  all  times  Np  =  0.  This  shortcoming  is 
corrected  by  adding  a  spontaneous  emission,  noise,  term  to  equation  1.  However,  it  is 
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simpler  to  deal  with  circulating  intensities  within  the  cavity  rather  than  photon  numbers. 
The  two  way  average  intensity  is  related  to  the  number  of  photons  by 


(2) 


where  V  is  the  mode  volume.  Rewriting  equation  1  in  terms  of  intensity  and  adding  the 
spontaneous  emission  term  results  in 


dt  d  VK  ’  /,h>  v 


(3) 


where  the  last  term  is  the  spontaneous  emission.  Gain  is  a  function  of  the  population 
inversion  and  was  given  as  equation  2.57,  and  is  repeated  here  in  terms  of  the  CO 
rovibrational  inversions  for  both  P  and  R  branch  lasing  transitions. 


rM  =  <r„  (v.  Nt?-x-J±x) 


(4) 


This  approach  implies  that  the  rate  equations  for  the  population  inversions  must  be 
developed  in  order  to  calculate  the  circulating  intensity  in  equation  3. 

Assuming  no  longitudinal  dependence,  the  energy  is  absorbed  uniformly  in  the 
pumped  volume  which  leads  to  the  following  source  terms  for  the  2nd  overtone 
population  [32], 


E  abs 

h  v 


J_ 

V 


(5) 


Where  Eabs  is  the  absorbed  energy  of  the  pump  beam,  and  V  is  the  mode  volume.  The 


time  dependent  source  terms  assuming  a  Gaussian  pump  pulse  is  described  by 

f  ' 


WUpper(t)  =  [Nvl- 


rp  V  n 


-41n(2)|^-‘0 


L 

p  j 


(6) 
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where  rp  is  the  pump  pulse’s  temporal  width,  and  t0  is  the  location  of  the  peak  in  time. 

Equation  6  is  the  only  source  term  for  the  2nd  overtone  manifold.  Each  rotational  level 
will  have  its  own  rate  equation  describing  the  population  as  a  function  of  time.  Each  rate 
equation  will  include  RR,  W,  VT,  stimulated  emission,  and  spontaneous  emission 
source  and  sink  terms  for  all  allowed  transitions.  Only  RR  will  be  included  in  the  model 
since  W  and  VT  take  place  on  time  scales  much  larger  than  the  lasing.  Only  the  specific 
rovibrational  pumped  will  have  pumped  source  terms,  which  drives  the  remainder  of  the 
system.  The  generic  rate  equation  for  the  pumped  level  and  non-pumped  level  are  given 
by: 

=wvr,r(t)-rA'’,J)-rA'’J)-4N(v,J)}+B(v,j)  (8) 

dt  Pump 

^,J)  =  Yp{v  +  hJ  -l)+  rR(v  +  l,J  +  \)-yP{v,J)-yR(vj)  ^ 

-A[N(v,j)]+B{v,j) 

=  yP  (v  + 1,  J  - 1)  +  yR  (v  +  i,  J  + l)  +  B{v,  J )  -  fVLower  (t)  (10) 

ut  Ground 

where  the  B(v,  J )  terms  represent  rotational  relaxation  into  and  out  of  a  rotational  level, 
and  J[jV(v,/)]  represents  spontaneous  emission.  Equation  8  is  the  rate  equation  for  the 
specifically  pump  rovibrational  level,  equation  9  is  for  all  other  rovibrational  levels,  and 
equation  10  is  represents  the  rovibrational  level  pumped  from. 

The  model  was  run  under  a  variety  of  conditions  such  as:  including  all  kinetics  or 
only  a  select  few  kinetic  terms.  It  was  observed  that  the  numerical  model  produced 
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similar  results  with  and  without  W  and  VT  kinetics  included.  Thus,  the  model  was 
typically  ran  without  VV  and  VT  to  decrease  the  computation  time.  Figure  C.l  has  an 
example  of  data  collected  with  various  effects  excluded. 


CO@95Torrw/  CM_Y  RR  WNET1CS 

15mJ  Pulse 
PIWP  QA:  0.30  cm 
CELL  L  =44" 

Rssonator  L  =  70" 


CUput  htensity  (P7  V=1) 
Cutput  htensity  (P9  V=1) 
CUput  htensity  (P8  V=2) 
CutpU  htensity  (R5  V=1) 
CUpU  htensity  (R7  V=1) 
CutpU  htensity  (R6  V=2) 
FLrrp  Rise 


Figure  C.  1.  The  negligible  effects  of  W  and  VT  kinetics. 


The  model  lacked  a  few  pertinent  parameters  such  as  a  longitudinal  and  radial 
spatial  dependence.  The  gain  was  considered  spatially  uniform  in  the  model  which  was 
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not  the  case  experimentally.  The  model  also  lacked  an  ASE  or  super  radiance  term 
which  was  a  dominant  factor  in  the  experimental  data.  The  final  improvement  would  be 
to  include  a  model  of  the  radiation  in  the  cavity  so  the  transverse  and  axial  modes  could 
be  accounted  for,  and  mode  locking  effects  could  be  modeled.  A  typical  model  output 
plot  is  in  Figure  1.  It  does  not  show  any  of  the  typical  temporal  variations  of  the 
experimentally  collected  data,  since  mode  beating  effects  are  not  included  in  the  model. 


Figure  C.2.  Model  output  of  an  v(0,3)R(6)  pump  case  closely  matching  experimental 
run  conditions. 

The  model  does  not  accurately  predict  the  temporal  output  of  the  CO  laser  due  to 
the  missing  attributes  mentioned  above.  However,  the  model  reasonably  predicts  overall 
trends  in  the  CO  laser  efficiency  versus  pressure,  output  coupler  reflectivity,  and  path 
length  effects. 
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